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Note on a Class of Exact Solutions 
in Magneto-hydrodynamics 


CeCeLiw 


1. Introduction 


The purpose of this note is to describe a class of exact solutions of the magneto- 
hydrodynamic equations which are applicable to various flow configurations. 
Some examples of these solutions have been worked out in some detail by PETER 
GOTTLIEB; his solutions and some others will be reported in future communi- 
cations. 

The class of solutions considered here is of the “layer-type’’, including many 
classical exact solutions of the Navier-Stokes equations as special cases; e.g., 
the flow near a stagnation point in the two-dimensional and the axially symmetri- 
cal cases, and the flow of a viscous fluid over a rotating disc. It represents a 
further generalization of a class of solutions discussed by the present writer for 
the steady or unsteady motion of an ordinary viscous fluid. 


2. Basic equations 
For a uniform electrically conducting incompressible fluid, the equations of 
magneto-hydrodynamics are as follows*: 


OU; 00; ic 0H; oll OUR 
2 ee peat yo Redes Ble Se, | CU hese 
(2745) = =o ame oT ee vA v,;, ax, ; 
oH; 0H; 00; oH 
2 P oY . U Fie +, — H. Boe siete 
2) of iby, i Ox, BAC 2 aia 


In these equations, x; (1 =0, 1, 2) are the space coordinates, ¢ is the time, v; are 
the components of velocity, H; are the components of the magnetic field. The 
properties of the fluid are characterized by the density @, the kinematic viscosity », 
the magnetic permeability ~ and the electric conductivity o, in terms of which 
we define the magnetic diffusivity 


(2.3) n = (4rpeo) 7. 
The quantity JT is essentially the total pressure and is defined 
(2.4) oll=p+yH*/8x, 


where # is the ordinary pressure of the fluid. The summation convention is used 
throughout this paper. 


* Cf. COWLING (1957), P. 93 
Arch. Rational Mech. Anal., Vol. 1 27 
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We seek a class of solutions in which the x-axis occupies a special position. 
We allow the velocity field, the magnetic field, and the pressure gradient to 
depend on the coordinate %» and the time ¢ in a fairly general manner, but only 
linearly on the two other coordinates x, («= 1, 2). To be more specific, we assume 


(2.5) Up = Uy (Xo, #) 
: Vy = Ug (%q, t) + Uap (Xo, 2) XB » 
(2.6) Hy = hg (%o, #) 
: Hy = hy (%, t) + hap (Xo, t) Xa, 
and 
oll a 
ay = Bo (%> ?) 
(2.7) ae 
ae G(X, t) + Dz,p (Xo; t) Xp 
Since we must have 
ra] ell \ ra) ell 
(2.8) 0%, ( aa = O%o i) ; 


the functions @,(% 9, !) and @,,(%9, t) must actually be independent of x). The 
function J/7 then takes on the form 


(2.9) IT = & ®z8(t) %q_ Xp + Ga (t) xe + D(x, 2), 


3. Equations in two independent variables 


The particular choice of the form of the solutions is made so that the other 
terms in the basic equations (2.1) and (2.2) are also linear in the variables x, 


(x=1, 2). Consider, for example, v; gis By using (2.5) and (2.6), it is easily 


0x; 
seen that 
0H; hy 
B-1) ere ys 
and 
0H ah ah 
(3.2) 2, a aa (2 ath hay) Bid (v4 Shas +2, a hay) Xp. 


Thus, these basic equations can be satisfied by comparing coefficients of x, and 
%, and terms independent of them. 


For convenience of reference, we list below the formulae useful for such reduc- 
tions. Each of the vectors v, and H; is of the form 
(3.3) {Aj} = {a@p, a, + dag Xp}, 


where dy, 4,, 4,, are functions of x) and ¢. The vector 0A,/0¢ is then of the same 
form. In fact 


(0A;\ _ f{ 0a, Ga, | Gay 
oe { a hay eg oe xp. 
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Also, the Laplacian of A; is of the form 


i DNGhy Gres O* Aap 
(3.5) {A A;} ={ ax ’ ox? =F Ox2 xh. 
Furthermore, if B; is another vector of the same form as A;, then 
ee ce 
is also of the same form; indeed, 
0a Oay 
Caer Pra Cy = dy ay my Oe 
6.7 i 
Cap = Oo oy, Fay Oyp 


With the help of these formulae, it is not difficult to write down the equations 
for the functions of x» and ¢ defined in (2.5) and (2.6). 


The equation of continuity in (2.1) is satisfied if 


OU as 
(I) Ox» + Ugg = 0. 


Similarly, the condition for the absence of free magnetic charges is satisfied if 


(1) Sat + hap = 0. 

The equations of motion (2.1) lead to 

gata) ay Ste Be +l oS), 

(IIIb) oe + Ug oe | Way Uae, oe on Dy vies ( Le | hay hy), 
(Io) “gE gt tay tye — 9 Gage = — Ban + Ghee (ha Gee + ay ya) 


Similarly, the equations (2.2) for the change of magnetic field lead to 


Oho |, oho Ohy 1, OUy 
yg) Ob yO Sxl nagar era? 

ah oh, Phy 1 Oly 
(IVb) aa + Uy ae tao axe 0 Bx, + Ury My, 

Oh, B Oh. Phe g nie OUg B 
(IV c) oF + Uo axe eee) Hyp — a8 0 Bx, + Un My g- 


4. Consistency of the system 


It is necessary to examine the above system of equations to see that there is 
indeed the right number of equations for the unknowns specified in (2.5), (2.6) 
and (2.7). Clearly, there are seven unknowns in each of (2.5) and (2.6) and only 
one unknown, 7.¢., @)(%o, ¢), in (2.7), since @, and @,g are functions of ¢ alone 
and are therefore expected to be specified by boundary conditions at specified 
values of xy (e.g. at infinity). Thus, there are /7fteen unknowns altogether. 


Duyfe: 
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On the other hand, we have one equation for each of the types (I) and (II) 
and seven equations for each of the types (III) and (IV), making up a total of 
sixteen equations. Thus, one of them must be redundant if we are not going to 
have the difficulty of too many equations. Indeed, a closer examination shows 
that the equations (I), (II), (IVa) and (IVc) are not all independent. 


If we contract (IVc), we obtain 


hep = Cups 


oh oh 
BB BB ay eB: 
at bt a bg is 1] axe “aso 0X» + Upy hyp, 
or 
, ON Bp Chee Phe p I OUg B 
a) af ty! oes ax 0 Bx, 


In this form, it is clear that, by using (I) and (II), «gg and hg, can be eliminated 
to give an equation in fy and wu). The resultant equation must then be consistent 
with (IVa). This fact can be easily verified by noting that the x)-derivative of 
(IVa), when added onto (IVc’), leads to an identity by virtue of (I) and (IJ). 


5. Concluding Remarks 


We have thus found a large class of exact solutions of the hydromagnetic 
equations, involving fifteen unspecified functions of the variables x) and ¢. If 
one restricts oneself to the steady case, the equations above reduce to ordinary 
differential equations which can be integrated (numerically, if necessary) under 
appropriate boundary conditions. Symmetry conditions also help to simplify 
the equations in various cases. 

These simplifications are all present in the examples worked out by PETER 
GOTTLIEB. He considered the flow near a stagnation point both in the two- 
dimensional case and in the case of axial symmetry with a magnetic field parallel 
to a solid boundary. These examples serve to illustrate the behavior of the 
magnetic field as it is carried along a solid surface due to convection currents in 
the core of the earth. It is clear, from a consideration of the path of the fluid 
particles, that there can be a considerable stretching of the magnetic lines. In 
the ideal case of zero viscosity and zero magnetic diffusivity, the field may be 
expected to become infinite as the boundary is approached. It is to be expected 
that this tendency to increase will be counteracted by both the hydrodynamic 
and the magnetic diffusive effects. The resultant picture in the case of steady 
motion is worked out in a paper to follow. Other interesting problems may be 
developed within the present framework. For example, one might consider an 
infinite disc oscillating about its axis (taken to be the x-axis) in the presence of 
a magnetic field in that direction. Some hydromagnetic damping is certainly 
to be expected. The steady rotation of a disc would be another possible problem, 
extending that worked out by von KARMAN & CocHRAN years ago* in the purely 
hydrodynamic case. Indeed, this may be combined with the problems studied 
by Gorr.ies to describe the combined effect of earth rotation and convection 
near a stagnation point in the problem of geomagnetism. Finally, one might 
consider the problem of a uniform magnetic field perpendicular to the Xg-aXiS 


* See GOLDSTEIN (1938), p. 110. 
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and rotating about it (perhaps in an unsteady manner), and attempt to find out 
whether this would indeed provide an effective “‘magnetic wall’ to hold an ap- 
proaching fluid. It is clear from elementary considerations of Eq. (2.1) and (2.2) 
that a magnetic field can act in the manner of a static pressure to hold a fluid 
in balance, but the added effect of rotation may have the advantage of making 
the equilibrium more stable. These and other problems are subjects for future 
studies. 
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Work Functions in Hypo-elasticity 


B. BERNSTEIN & J. L. ERICKSEN 


1. Introduction 


According to linear elasticity theory, the work per unit volume done by the 
stress in deforming a material is uniquely determined by the initial and final 
states of stress, regardless of what intermediate stresses it may experience. The 
theory of hypo-elasticity proposed by TRUESDELL [1J, 2, 3] constitutes one gener- 
alization of linear elasticity theory. As will be made clear later, the work density 
for hypo-elastic materials generally depends on the entire stress history. In 
other words, it is a functional rather than a function of the stress. The main 
purpose of this paper is to characterize two exceptional types of hypo-elastic 
materials for which the work density depends only on the end states of stress. 
Tuomas [4, 5] has given examples of materials similar to those discussed here. 
The materials of type II discussed below are included in the class of materials 
characterized in [6]. The work of CapRIOLI [7] suggested to us some of the results 
presented, though his analysis is not directly applicable to hypo-elasticity. 

We begin by establishing theorems on functionals which are useful for our 
purpose. The remainder of the paper is devoted to the application and inter- 
pretation of these theorems in hypo-elasticity theory. 


2. Functionals considered 
We study functionals of the form 
aye fie (yon) 22 ay 
roae.d= f ive] 4 vm tek PRL (2.1) 
0 
where Greek indices run from 1 to m, the summation convention is used, F, 
and G, are specified functions of y’, and y* = y*(z) represents an arbitrary curve. 
Except when otherwise noted, it is assumed that F, and G, are continuously 
differentiable functions of y’ in some region St and that the curves considered 
lie in ® and are continuous and piecewise continuously differentiable functions 
of 7 in the t-intervals considered. Results not clearly labelled as local are global, 
the conclusions holding throughout any subregion of 9 in which the hypotheses 
hold. For ease in writing, we abbreviate (2.1) by the expression 


Tt 
S Gap? (n) 
0 


F[j*(x),]= f Rape 


when the curve is given by y*= f*(z). 


(2.2) 
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Lemma 1. If the curves y*=f*(t) and y*=g*(t) are related by 


tm Pla—) =e(0), (2.3) 
in the t-interval <0, t), then 
ty 
- = f Gygax) 
F[f*(t),4,])=—e ° Tie), t\\- (2.4) 
Proof. In the left side of (2.4), make the changes of variable 
T=h—t, H=4—n. (2.5) 
obtaining by (2.2) 
. 0 ie apa) 
FIP @,4l= | Rap @e 
ty 
— f Ggahm 7 _ f6,ar@ 
=—¢ f Rap@e 
0 
— f Gah @) c 
=—2 9” * FIP, A). 


Using (2.3) and (2.5), we then obtain (2.4). 
Lemma 2. Let y*=a*C ht represent any fixed point and let t, denote any fixed 
t-value. If F[f*(t), 4,]|20 for every curve y*=f*(t) such that 
LO) ea IA Gh) "a, (2.6) 
then F [f*(t), 4] =0 for every such curve. 


Proof. Given any curve y*=/*(t) which satisfies (2.6), we can construct 
the curve 


which also satisfies (2.6). By hypothesis, 
F[g*(t),4)20, FUP (x), 4) 20. 
From Lemma 1, this is impossible unless F [/*(t), 4] =0. 
CAPRIOLI [7] established a lemma which is equivalent to our Lemma 2 when 
G,= 0. 
Lemma 3. Let y*= a" represent any fixed point in KR. A necessary and suffi- 
cient condition that, for each choice of t*=0, 


F[jf*(),t*]=0 whenever fP(0) =P (*)=a& (2a) 
is that there exist a function y=y(y"*), 
y(2*) =0, (2.8) 
such that, for every t 20, 
t 
fGgaf(x) 
F[f*(z),t]=ylP(@)]e (2.9) 


whenever 
f{O =a (2.10) 
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Proof. From (2.8), (2.9) and (2.10), we have 


f Gain J GyasAee 
F [f* (zt), *) =y[P (4)] & = y(a*) =0 


whenever /*(0) =/*(¢*) =0, so these equations imply (2.7). Suppose now that 
(2.7) holds. Select any point y*=6*¢ }i and any two positive numbers ¢, and fy 
and consider any two curves 
y=fx(t), OStSiy (N =4,2) 
such that 
fy (0) = 4", — fy (tw) = 8. (2.11) 


When b*=a%, we may take 4:=0. The curve y*=/*(t), where 
P(t) =A), O84, 


(2.12) 
=f4+4—1) AStS4rh, 
satisfies 
PO=P)=4, M=4+8h. 
From (2.7), (2.2) and (2.12), 
ttt, B 
fo df? (n) 
0=Fif*(2),4 +4] = J F,df*(t)e * 
t ttt, 
: f Ga en so Gp a/P(n) 
= [ Rap(e ° wel Fd f*(t (2.13) 
0 


p2fP(n) aS fc, af?(n) 


= Fifi(0), pega f Eafe 


Introducing the changes of variable t =4,+?,—1, 7 =4,+t,—y and using (2.12), 
we obtain 


t+, G, af” ( (7 . af? (7) 
f Fd f* (zt) ee = [Rae ft U7] 


t 


=—e 0 


fic,ai (7) {svat 
fr 2473) 


eee FUfg(t), te], 
which, with (2.13), yields 


ty 2 
— f Gap — f Gath 
ens Layee | eeroene BL) (2.14) 
This must hold for every pair of curves satisfying (2.11), where 0%, #, and ¢, are 


assigned arbitrarily. Thus a function y(y%) is well defined by 


= J Gye on 
y(y*)=e ° lb ea ca eal (2.15) 


—————— 
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where /*(t) represents any curve passing through the point y*such that /* (0) =a’, ¢’ 
being a value of t taken on when the curve passes through y*. For any such curve, 


(a) =y[P (0)] =F LP (x), 0] =0, 


which establishes (2.8). From (2.15), (2.9) follows immediately when (2.10) holds. 


Lemma 4. In order that there exist a function y(y*) and a point y*=a* such 
that 


t 
f GgafP(r) 


FIP (t), =v id) e (2.16) 


whenever {*(0) =a% it is necessary and sufficient that there exist a solution y(y*) 
of the system 


oy 
Bye Ga (2.17) 


such that, at some point y*=a*, 
Wri hae) eee (2.18) 


Proof. Suppose (2.16) holds. For any continuously differentiable curve 
y*= f*(t) such that /*(0) =a”, y[f*(é] is then a differentiable function of ¢ and 


t 
ee, ae ajx) {Sah 
Er oe a Oineyere (2.19) 


From (2.2), we also have 


d 6 og & 
qe lt (t), 7] = Lia e : (2.20) 
so that 
dy d{* d{* 
le nd gpa a ae 224) 


ahah 3 d f{* 
Now /*(0) =a* imposes no restriction on the values of = for f=0 or on the 


values of /*(¢) for £>0. It thus follows from (2.21) that the partial derivatives 


pee exist and satisfy (2.17). That (2.18) must hold follows from (2.16) and the 


fact that F [/*(z), 0] =0 for all curves. 

Suppose conversely that (2.17) and (2.18) hold. Then, for all differentiable 
curves, (2.20) and (2.21) hold, from which (2.19) follows. Integrating (2.19) from 
t=O to ¢, then using (2.18), we obtain (2.16) whenever /* (0) =a. 


Lemma 5. In order that (2.17), considered as a system of differential equations 
in the unknown y, admit more than one solution, it 1s necessary and sufficient that 
there exist, at least locally, functions A(y*) and O(y%) such that 


G,Son/oy", (2.22) 
F, =e“ 00/dy*. (2.23) 


Proof. It is easily verified that, if (2.22) and (2.23) hold, (2.17) admits the 
solutions y =e “(9+ A), where A is an arbitrary constant. Conversely, if (2.17) 
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admits two distinct solutions y, and y., the difference ¥1—72=A satisfies 


0A 
Gy = Gis, 
from which (2.22) holds with A=—InA. From (2.17) and (2.22), 
OV oA one A 


LL Bat Vege ages 
so (2.23) holds with O =», e*. 


3. Functionals of type I 

We pay special attention to two types of functionals, the first being charac- 
terized as follows: 

Definition 1. A functional of the form (2.2) is of typeI provided there exists 
a point y*=a%ER such that F[f*(z), t]20 for any t20 and for every curve 
y*= {* (rt) such that {* (0) =a*. 

Theorem 1. A functional is of typel if and only if there exists a point 
y*=a* ER and a function y(y*) such that 


wly")=0, yla)=0, (3.1) 
and such that, for any t=0, 


J Gpai®en 
Fif(t), tJ=rP Ole (3.2) 
whenever f*(0) =a. 3.3) 
Proof. From (3.1), (3.2) and (3.3), it follows immediately that F [/*(t), 4] 20 


whenever (3.3) holds, so these conditions are sufficient for a functional to be of 
type I. 

For any functional of type I, the hypotheses of Lemma 2 are satisfied for 
arbitrary 4,20. From Lemma2 and Lemma}, there then exists a function 
such that (3.2) holds whenever (3.3) does. That y must be non-negative follows 
from (3.2) and Definition 1. 


Corollary 1. A functional is of typel if and only if there exists a point 
y*—= a* and a function y such that 


Byte — 7 Ga tf (3.4) 


and such that us f 
y(y*)20, ya) =0. (3.5) 


Proof. This result follows immediately from Lemma 4 and Theorem 1. 


Corollary 2. For any functional of type I, 
F,(a°) =0. (3.6) 


Proof. From (3.5), y takes on zero as a minimum value at y*=a* so that 
dy/0y*=y =0 when y*=a*. From (3.4), we then have (3.6). 

Using Corollary 1, it is easy to construct functionals of type I: we can choose 
the functions G, arbitrarily, select any non-negative function y which vanishes 
at some point, then use (3.4) to calculate F.. 
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4. Functionals of type II 

The second type of functionals to be considered is described as follows: 

Definition 2. A functional of the form (2.2) is a functional of type Il provided 
that, for each t,=0, F [f*(t), t,] 0 for every curve y*= f*(t) such that f* (0) = f* (4). 

Theorem 2. I} the trivial case F,=0 is excluded, the following three statements 
are equivalent: 

a) The functional given by (2.2) is of type II. 

b) For each t,=0, F[f*(t), 4] =0 for every curve y*=]*(t) such that f*(0) = 
(4). 

c) There exist functions O(y*) and A(y*) such that, for any t=0 and for every 
curve y*= f*(z), 


[GaP () = Al] — AIO), (4.1) 
F [f* (x), f] = {O [f* ()] — O [f* (0) J} e407, (4.2) 


Proof. That a) and b) are equivalent follows immediately from Lemma 2 
and Definition 2. From (4.2), it is clear that c) implies b). It remains to show 
that b) implies c). When b) holds, it follows from Lemma 3 that, for each 
choice of a%, there exists a function y(y*) such that, for any ‘20, and any curve 


y= f(t), t 


ad Gp f(x) 
FiP(t), J=yPO] hile 0 
whenever /*(0) =a*. Equivalently, there exists a function y of pairs of values 
of y* such that, for any ‘0 and any curve y*=/*(t), 
t 


f Gp a P(x) 
FF ().=71F0.P Ole? , 
or 
= “Gga? 
PO) eer? F [f*(z), ¢] ay 
: Gp afP(n) é 
= f Rape , ¢{f"(0), f”(0)] = 0. 


Select any three points y*=a*, 6% and c* in ® and any curve y*=/"(t) such 
that, for some t values ¢, and 4,240, /*(0)=a%, /*(4) =0* and /*(é,) =c’. 
From (4.3), we then have 


t 
fGgah(n) 


ec) = f Edf (ne 


J dfP(n) fe,a%m fe,acm 
Son ‘J Rare a” +f Rafe * 
ty 
'GpdfP(n) tot i) GpafP(r) 
oi B Tae aX 0 y+ i dr ( aos wes: A 
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where T=1t—4,, 7 = —t,. From (4.3), the last term has the value VAUBRECS \- 
-We thus have t, 


Sf GgafP(x) 
x, (a, c”) = ef 7 (a%, 6”) + 4(0%,¢"), (4.4) 
ay PO=% Fah, Fi) (4.5) 
Since /*(0) is arbitrary in (4.3), 
y(y*,y’) =0~ for all y*. (4.6) 


When EF, does not vanish identically, there always exist pairs of points for which 
y~+0, eg., two distinct points on a curve satisfying dy*/dt =F, one being a 
point at which F,--0. Let a* and b* be two such points and consider the case 
where c*= b* but #, and ¢, may differ. From (4.4), (4.5) and (4.6), we then have 


ty 
SG,af()=0 whenever PU)—=-) =o (4.7) 
from which (4.1) follows immediately. From (4.1), (4.4) and (4.5), we then have, 
for allva? 0 andes 
x (a*, c”) ef (os) eas 4 (a%, b”) ef! (64) x (b%, c”) ef! (cM) (4.8) 
Setting * a7==0,.0'(C) =» (0567) & (cP), we obtain from (4.8) 
eA (") 7 (b%, c”) = O(c") — O68"), (4.9) 
for all 6% and c*. From (4.3), (4.1) and (4.9) 
t 
q Gy 4f%(n) 
F[f*(t),t] =e xP (0), PO] 
= 6 ee OU at 
= 6 AMO {0 [f° ()] — OL ()]}, 
which is (4.2). Hence b) implies c). 


The functions @ and A are not uniquely determined by (4.1) and (4.2). As 
is easily verified, we have 


Leo oe 4.10 
e410") — OW] =e *[0'(y4) —O"E*)], ah 
for all y* and 2% if and only if there exist constants A and B such that 

A=A'+A, O=c4(O'+ 8B). (4.41) 


Theorem 3. The functional given by (2.2) is of type ll if and only if there 
exist functions A(y*) and O(y%) such that 


G, = dAlay*, (4.12) 
E= 2 S00loy™ (4.13) 


* If a*=0 is not in §t, we may select any other fixed point which is. 
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Proof. From Theorem 2, it suffices to show that (4.1) and (4.2) hold if 
and only if (4.12) and (4.13) hold. Clearly (4.1) holds if and only if (4.12) holds. 
Differentiation of (4.2) with respect to ¢ gives, for arbitrary continuously dif- 
ferentiable curves, i 

dpe [GaP oe 
fo, @ ee a7 Ol). (4.14) 


Using (4.1), we may write (4.14) in the form 


Oe B d 
Bone. = aero, (4.15) 


from which (4.13) follows. Conversely, if (4.13) holds, (4.15) follows for any 
continuously differentiable curve. Multiplying (4.15) by e~4%) and using 
(4.1), we obtain (4.14). Integration of (4.14) yields (4.2). 


From Theorem 3, Lemma 5 and Corollary 1, a functional of type II is 
similar to a functional of type I in that, for either type, there exists at least one 
function y satisfying (3.4). For a functional of type II, it is not necessarily true 
that there will exist a solution of (3.4) satisfying (3.5). 


Theorem 4. For a functional to be of typeI and of type Ul, tt is necessary 
and sufficient that there exist functions A(y*) and O(y*) satisfying (4.12) and (4.13) 
such that © takes on an absolute minimum at some point y* =a* CH. 


Proof. From Theorem 3, a functional is of type II if and only if there 
exist functions A and @ such that (4.12) and (4.13) hold. From Corollary 1, 
a functional is of type I if and only if (3.4) and (3.5) hold for some point y* = a*. 
A functional is thus of both types if and only if there exist functions y, A and O 


such that 
000 dy oy 0A 
Ee =—€ : i, yo Fr ay = YW Gy One iz i ay x (4.16) 


with y (a) =0, y (y*) = (@) (4.17) 


for some point y*=a*. For (4.16) to hold, it is necessary and sufficient that 
(4.12) and (4.13) hold and that 


y=e“(@+A), A=const. (4.18) 
For (4.17) and (4.18) to be consistent, it is necessary and sufficient that 
Oly") 20(@"), A= —O(a"), 
which completes the proof. 
5. Hypo-elasticity equations 


Following Nott [9] and Tuomas [10, 11], we write the constitutive equations 
for hypo-elastic materials in the form 


jn v7 ; coe 
t= ] — biz Win — bp Wir = Aijnm tam (5.1) 


using Cartesian tensor notation, where ¢;;=¢;; represents the stress tensor, ¢ the 
time, d/dt the material derivative, w;, and d;, are given in terms of the velocity 
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components v,; by 2@;,=%j,4—Ux,7) 24x = Yj, b+, j and A;;,,, represents a tensor 
invariant or eanropie function of ¢,, with the symmetry A jj%m=Ajizm = Assman 
Dt,;/dt transforms as a second order tensor under arbitrary time-dependent 
orthogonal coordinate transformations and reduces to 0¢,,/8¢ at a point P in a 
reference frame such that v;=w;,=0 at P. We also use the equation of con- 
tinuity in the equivalent hai 


d 
+ 0dn=9, dV =o dh, (5.2) 


o and V being the density and volume occupied at time ¢ by a set of particles 


initially occupying the volume VY, and having density gp. We require that (5.1) 


may be inverted to give a 
km 


d;;= Bijzm Ge (5.3) 


in some region 9 of stress values, B;;,,, being a tensor invariant of ¢,, with the 
same symmetry as A;;,,,- Here B;;,,, Which is the inverse of 4;;,,,, considered 
as a 6X6 matrix, 


Aicnm Deans <a Bes has Are f- z (6;, 0; . 1% 0; 0; ,) , (5.4) 


exists if and only if D#,;/dt = 0 implies d,,,=0. Since B;;,,, has A;;,, aS an inverse, 
In KR, B; jm Or, equivalently, A;;;,, 1S assumed to be a continuously differentiable 
function of ¢,,. We henceforth exclude stress values not in #. Continuity 
requirements for the stress and motion should be evident from the context. 

We find it convenient to introduce as independent variables the time and 


material coordinates X;, such that the motions considered are represented by 
relations of the form 


Hi = %j(Xy,2), Xy= Xp, (%j,8),° 4; (Xp, 0) = Xj. (5.6) 
Then a particle is represented by fixed values of the X,’s and 


ad B. Qi Bo ax; 
de OR Leow Vi aes 


It should be noted that, under the above assumptions, (5.1) and (5.2) impose 
no conditions on how the stress depends on the time at a single particle: we can 
always choose a stress depending only on the time which coincides with any 
given stress at a given particle at all times, arbitrarily assign functions W;; = — Wj; 
of the time only and a positive constant 9, use (5.3) to calculate i (5:2)"t 

5: 


calculate @, and take v; = (w;;-+d,,) x;. With this assignment, (5.1) and (5.2) are 
satisfied. 


6. Work and density functionals 


The work W(t) done by the stress in deforming a material volume V (é) during 
the time interval <0, t> is defined by 


Sf tid (x) dt. (6.1) 
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It is noted in [6] that dt;; 
whenever (5.1) holds, where 
Fj = Bu mij m= Fi; (6.3) 
From (5.2), t 
0 f denae 
EV ego avy, (6.4) 


where the integral is a partial integration with respect to ¢ holding X, constant. 
From (5.3), 


DE 
dy =Gyj 7, (6.5) 


where 
G,;= Brig = Gi ; (6.6) 


is an isotropic function of ¢,,. As such, it is expressible* in the form 
Making this substitution in ae 6 and using (5.1), we see that 


dt;; 
di, = Gia, (6.7) 
the terms involving w,; dropping out in the summation. From (6.4) and (6.5) 
t 
f Giz ati; (2) 
Sa (6.8) 


which expresses the density at X, as a function of its initial value 9, and a 
functional of the stress at X,. We also have from (6.1), (6.2), (6.4) and (6.8) 


ef Jia pp aNydt = | w(X,,0) d%, (6.9) 
Vy 
where 
0 [Gimatin'n ) 
wed) = ftv’ ips [i yatt0 . (6.10) 


We thus have the work per unit initial volume at X, expressed as a functional 
of the stress at X,. The tensors F;;, G;; and ¢;,, being symmetric, have six 
independent components. With the Gace 


ein mar act ae yp, y®) = (aes too, 33,3 (419+ f01), & (tes + f30)) 3 (43+ ts 1)) ) 
(fi, FF, Ea) (Ain, Foe, 3, 2F2, 243,243), (6.11) 
(G,, Gz, Gg, Gq, Gs, Ge) = (Gir, Goo, Gag, 2Gi2,2G23,2G;3), 
we can write the functional in the form (2.2). Henceforth 6£/0t;;, where & is 
any function of the stress, shall mean 


gE _ OF l1,, oz 
Gy; me Oh; 2 (nn 4 


(6.12) 


* See, e.g. [8, Appendix I], setting aj) =14;;, a}? =O. 
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where the nine quantities ¢,; are regarded as independent in performing the 
differentiation. With (6.11) and this convention we have, for example, the 
equivalence of the equations 

oA oA 


Ga pao 5 Oe rate 


7. Materials of type I 


Any assignment of the functions A;;,,, subject to the conditions mentioned 
in Section 5 serves to define a hypo-elastic material. We call such a material 
a material of type 1 provided its work functional, given by (6.9), is a functional 
of type I. This amounts to saying that a hypo-elastic material is of type I pro- 
vided there exists a stress value 7}, such that the work W(i) done by the stress 
in deforming any material volume V(#) is non-negative whenever ¢;; =¢?; at t=0, 
this being true for all motions and stresses consistent with the constitutive 
equations. 


We make free use of (6.11) and (6.12). 

From Theorem 1, a material is of type I if and only if there exists a stress 
value ¢?, and a non-negative function y(f;;) such that y(#,)=0 and such that, 
for any ¢=0 and fixed X,, 


t 
if Crm Ft m() 


w(t) = y[t;()] & (7 1) 


whenever ¢;;(0) =¢,, w(f) being given by (6.10). From Corollary 1, a material 
is of type Lif and only if there exists a stress value #7, and a non-negative func- 


tion y (¢;;) such that 
a] 
é. s = —»G;;+ F,, y (,) =0, (7.2) 


and, from Corollary 2, 


Theorem 5. lor any material of typeI, t?,=0, 
W(t) =f 7 Lee, t)|}dV(t) when #,;(X,,0) =0 (7.4) 
and, regardless of whether t;;(X;,,0) =0, 
d 
t;;4;; =@ dt (y/e) - (7.5) 
Proof. For any material of typeI, we have from (6.3) and (7.3) 
Fi; (8 m) am Bu mij (tq) tem Evy 


which, from (5.5), implies that #,,=0. From (7.1) and (6.8), for fixed X,, 


km 


w (t) = rl; O19 when ¢;;(0) = 0. (7.6) 
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Integrating (7.6) over Y), then using (5.2), we obtain (7.4). From (6.2), (7.2), 
(6.7) and (5.2), 


dt; ; Cy at; dt;, 

bij dij = Feige wen Tce Cr ap 
dy dy _ y do 
dt tue = Gi — ¢ dé’ 


which establishes (7.5). 


Thus y represents the work per unit present volume done by the stress in 
deforming the material from an initially unstressed state. 


Theorem 6. A hypo-elastic material is of typeI if there exists a function y (t;;) 
such that, whenever (5.1) and (5.2) hold, 
t54;;=oq Vie), yvh)20, yv(0)=0. ~ (7.7) 


Proof. Under these hypotheses, we have, using (6.10) 


w) =f td; tdr =f ogc tyledr=“y [4,4] =0, (78) 


whenever (5.1) and (5.2) hold and #;;(0)=0, X, being fixed. Since (5.1) and 
(5.2) do not restrict the values ¢;,;(t) for fixed X,, the functional w(é) is of type I 
by Definition 1, (7.7) and (7.8). 


Theorem 7. A necessary and sufficient condition that a hypo-elastic material 
of type I have the property that there exist, in the neighborhood of each stress value, 
stress values which cannot be attained by deforming the material subject to the 
condition 


is that there exist, at least locally, functions D(t;;) and V(t;;) such that 


ow 
E,=@ Bh (7.9) 
1 ow oy 


Proof. [t is shown in [6] that a hypo-elastic material has the indicated 
inaccessibility property if and only if there exist, at least locally, functions ® 
and Y of the stress such that (7.9) holds. For a material of type I, (7.9) and 
(7.10) are equivalent since (7.2) holds. 

When there exist functions ® and W satisfying (7.9), W is called a hypo- 
elastic potential. 

8. Materials of type II 

A hypo-elastic material is a material of type II provided the work functional 

w(X,, t) is a functional of type II. This is tantamount to requiring that, for all 


material volumes and for all stresses and motions consistent with the constitutive 
Arch. Rational Mech. Anal., Vol. 1 28 
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equations, the work W(d) done by the stresses in deforming the material be non- 
negative whenever the stress path is closed, 


ti)(Xp,t) = t,;(Xp, 0)- 
From Theorem 2, a material is of type II if and only if there exist functions 
A(t;;) and O(¢;;) such that, for any ‘20 and any stress path 

FG, dtj(0) = Alby (9) — Afb (0), (8.1) 
w (t) = {O[t,;(] — O[t,;(0) fe AM, (8.2) 

X, being held fixed. From Theorem 3, (8.1) and (8.2) can be replaced by 
G,; = 0Alat;;, (8.3) 
Ee é “dOl0t,,. (8.4) 


For any such material, F; =®dW/dt,,, where D=e"“*, Y=O, so a hypo- 
elastic potential exists. From Theorems 3 and 4, a material is of type I and of 
type II if and only if (8.3) and (8.4) hold and @ takes on an absolute minimum 
at ¢;; =0, in which case 


y (t;;) =e" 4"? (O@,,) — O)]. 
Also in this case, (7.9) and (7.10) hold with @=e*, W=@. 


Theorem 8. Choose any fixed stress value t?,; and any material of type II. 
There exists an admissible choice of A and O such that, whenever t;;(Xz, 0) =t;, 
we have 


W(t) Sh Ot; (Xz, 4) aN, (8.5) 
@= aoe“, (8.6) 

oe dO 
54s = Oo ae (8.7) 


Proof. As noted earlier, A and @ are subject to the transformations (4.11). 
By proper choice of « and 8, we can arrange that O(/?;) =A(é,) =0, where #, 
is specified. With this choice, we obtain (8.5) from (8.2) and (6.10) and (8.6) 
from (6.8) and (8.1) when ¢;;(X,,0)=¢;. From (6.2), (8.4) and (8.5) 


E dt;, eA 00 at; ; Pst O dO 


E So al ———--  —— a 4 = 
GEES Lie Wane Figo tues 


d 


vy 
which is (8.6). 

Thus @ represents the work per unit initial volume done by the stress in 
deforming the material from a specified state of stress when @ and A are suitably 
normalized. Under the same conditions, the work per unit present volume Q is 
given by 
D=0 oo, =e- 0, (8.8) 


both these work densities being functions of stress. 
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Steady Shear Flow of Non-Newtonian Fluids 


WILLIAM O, CRIMINALE, JR., J. L. ERICKSEN & G.L. FILBEY, JR: 


I. Introduction 


In recent years, considerable effort has been expended in attempting to gain 
an understanding of the behavior of non-Newtonian fluids in shear. MARKoviTZ[1] 
recently attempted to collect and organize experimental data and to compare 
it with predictions of existing theories. One difficulty is that conclusions which 
can be drawn from different types of experiments seem to disagree. 


Theoretical work on the steady flow of non-Newtonian fluids has brought 
to light a new phenomenon. According to the linear theory of viscous fluids, it 
is always possible for a fluid flowing through a cylindrical tube, to which it 
adheres, to undergo steady rectilinear motion, each particle moving with constant 
speed in a straight line parallel to the generators of the cylinder. Analyses made 
by ErIcKsEN [2] and by GREEN & RIVLIN [3] show that, for many ideal fluids 
described by the REINER [4]-RIVLIN [5] theory, this simple type of motion is 
possible only for very special shapes of tubes and that, in general, it is replaced 
by a flow consisting of such a rectilinear motion combined with a secondary 
flow in cross-sectional planes. Using the more general theory of fluids proposed 
by Riviin & ERICKSEN [6], LancLots [7] has obtained solutions involving 
secondary flows in tubes and other types of boundary geometries. As has been 
discussed by RIvLIN [8], similar secondary flows have been observed in the mean 
flow of fluids through non-circular tubes when the motion is turbulent. To our 
knowledge, it is not known whether secondary flows occur in the steady flow 
of any real fluid through tubes. | 


Our analysis indicates that, if the conclusions which Roserts [9] drew from 
his experiments are correct, the fluids which he observed should be capable of 
undergoing rectilinear motion through tubes of any shape. Markovitz [1] draws 
self-consistent conclusions from other experimental data which contradict some 
of RoBerts’ conclusions. Calculations made in accordance with these conclusions 
suggest that secondary flows are to be expected in some shapes of tubes. We 
offer a partial proof that if secondary flows do not occur in elliptical tubes, they 
will not occur in tubes of any other shape. 


II. Equations of Motion 


For incompressible fluids, the constitutive equations considered are of the form 


is 2)T fA (2.1) 
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where f€ is the stress tensor, p is an arbitrary pressure, N is any positive integer, 
the A’s are symmetric tensors with rectangular Cartesian components given by 


(Ay);;= Cig Uy x 
ra) 
(Ayiiij= ae (Am)i; + (Am)iz.n% + (Am)in2s,; + (Amini > (eZ) 
M = —-.,.N — 1, 


where v, are the velocity components, and f denotes any symmetric isotropic 
function of the A’s, to be determined by experiment. These quantities must 
satisfy the incompressibility condition 


Unk = O, (2.3) 
and the equilibrium equations 
Ov; 
A Sree ee oA 


where @g is the density and @ the potential of the body force per unit volume, 
assumed conservative. 
From (2.2) and (2.3), 
trAe=tr Asin At =: (2.5) 
Whenever 
Ay 0 ior re 2 or FN = 2. (2.6) 


f reduces to a symmetric isotropic function of A, and A,. As such, it is expres- 
sible in the form [6], [10] 


f = %1 + 0 Ay + %, Ag+ a3 Ai + ot, A3 + «5 (A, A, + A, Aj) + 
+ a, (4? A, + A, At) + a, (43 A, + A, AB) + 0 (A? AB + AB AD), 
where 
ar=op(trA,, tr Aj, tr Aj, tr A,, tr A}, tr A}, trA,A,, 


‘ ; oe (2.8) 
tr Aj, A,,trA,A , tr Aj Ab), 


the «’s being polynomials whenever each component of f is a polynomial in the 
components of A, and A,. A representation of the form (2.7) with ag=0 always 
exists. With it, the remaining «’s may not be polynomials if f is. Absorbing a 
in the arbitrary pressure p and using (2.5), (2.7) and (2.8), we may write (2.1) 
in the form 


t= — pl + %A,+ %Ay+ 0g At + 04 AS + 015 (A, Ay + Ap Ay) + 


2 (2.9) 
+ o%¢(Aj A, + Ay Ai) + «,(A3 A, + Ay Ab) + 0% (AA + A3 A}), 
where 
Cr= ap (tr Aj, tr Ay tr Air Ass tr AA; tA; A, tr A, A3, tr Aj A}), oe 
Bat nc. 98. i 


Using these equations in essentially this form, R1vi1n [17] has obtained solutions 
for simple shearing motion, torsional motion, and combined Couette-Poiseuille 
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flow. The Reiner-Rivlin theory is obtained by setting N=1 in (2.1). For it, 
(2.9) holds with : 
hy = % (tr Aj, tr Aj), OO a, (tr Aj, tr Aj), 


(2.11) 
Ap = Ag = As = Ag = A — Oy QO: 


III. Laminar Shear Flows 


We call a motion a laminar shear flow provided that it is possible to choose 
rectangular Cartesian coordinates at any point P such that, at P 


0 |X 0 0: =0->-0 
A= |x 0 o|, A=10 2X 0], Ay=—Oo' for M2 oe 
Ona OO Of Ome O 
where 
X = 7 trA}. (3.2) 


For any laminar shear flow, we have from (2.9) and (2.10) 


t= — p1+,A,+ B,A, + B3(Ai — $A), (3.3) 


where the f’s are functions of X, given by 


Br=1t 275X + 4, X 
Bo = Ye 3 V3 + 2(Yat ttn ack 


Bs = Ys; 
Vr A) Op 2X, 0,4 AO ro O 2X Oar le 


(3.4) 


These relations hold in all coordinate systems. MARkovitz & WILLIAMSON [12] 
have used the constitutive equations in a form very similar to and equivalent 
to (3.3). RIvLIN’s general solutions are all laminar shear flows. For a Reiner- 
Rivlin fluid, it follows from (2.11) and (3.4) that 


Bo = 2B3- (3.5) 


For the more general fluids considered here, 8, , 8. and f, can be assigned arbitrarily 
as functions of X. We assume that £,-+0. It is easily shown that (3.5) holds if 
and only if the principal directions of ¢ and A, coincide. 

We are primarily concerned with the rectilinear motions 


1 = f(%2,%3), V2 = V3=0, (3.6) 
which satisfy (2.3). From (2.2) and (3.6), 
OE ke OC 0 
AiG 0. Oh Ag 0 12 -fels|>. Aw= Ole eM 
hfe mae) O fete ts 


By a rotation of coordinates in the x,-x, plane, we can always make f , 


f3=0 at a point, so that (3.1) holds. Thus these are laminar shear flows oa 
that (3.3) holds. From (3.2) and (3.7), 


X= fiefs Shale (3.8) 
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We have from (3.7) and (3.3) 


0 fe fs One X00 
t=— 1+, he Oh0 + 2B, 0) te hel s + Bs 0 O | ‘ (3.9) 
fs OF 0) hots ies OFEO 550 
Substituting these stresses in (2.4) we obtain 
P= (Bif,:),i> 


P, 3 2(Bsf ef:),a> 
P= 2(Bofstha),i 
D>? + ¥. 


(3.10) 


These equations are of the same form as those analyzed in [2], from which 
we may conclude that, for a solution to exist, there must exist a constant a 
and a function R(f) such that 


a= (Bifi),a (3.11) 
FF = 2B sh.) (3.12) 
P=prp=R(f) + f pedX rar. (3.13) 


For a given fluid, a and R(f) may be different for different flows. If 
2 ok Dyn eee = CONSUL, (3.14) 
every solution of (3.11) clearly satisfies (3.12) with R(f)=kaf. If 
Gre m/|/X, m= const., (3.15) 


it follows from STONE’s work [13] that (3.11) implies that the curves /=const. 
are straight lines if a=0, circles of radius 1/|ma| if a+-0. If neither (3.14) nor 
(3.15) holds, it is shown in [2] that, for each choice of a, there exist solutions 
of (3.11) which do not satisfy (3.12) for any choice of the function R(/). 

We define a tube to be a cylinder, S, with generators parallel to the %,-axis 
whose curve of intersection, @, with a plane +, const. is a simple closed curve. 
For a fluid to be capable of undergoing steady rectilinear motion through the 
interior of S and of adhering to S, there must exist at least one solution of (3.10) 
inside @ such that 

f= Orx0n! X64 (3.16) 


For a Newtonian fluid, 8, = const., 6,=/;=0, and there exists a unique solution 
for any given constant a, whenever @ is sufficiently smooth. In this case, (3.14) 
holds with k=0. For more general fluids, it is to be expected that, when a 
solution exists, f will be determined essentially uniquely by (3.11), and that the 
solution thus obtained will satisfy (3.12) for some choice of R(f). A solution 
may fail to exist because f, is such that (3.11) admits no solution satisfying 
(3.16). As is easily shown, this is the case for all choices of @ if (3.15) holds. 
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This is the case for so-called perfectly plastic solids, which are formally included 
in the theory considered. If this difficulty is not encountered, a solution may 
fail to exist because the solution(s) determined by (3.11) is(are) not consistent 
with (3.12). From Riviin’s analysis of Poiseuille flow [//], it would appear 
that this latter difficulty rarely, if ever, occurs when @ is a circle. From the 
remarks made above, it is to be expected that this difficulty will be encountered 
for many choices of @ if (3.14) does not hold for the fluid considered. The analyses 
made in [3] and [7] indicate that the assumed rectilinear motion is then replaced 
by a motion of the type 


0, =f (%2, %3), Y2—= Ya» Vo ee (3.17) 
y= (Xa, %), 
satisfying the adherence condition 
Ui, 0) Ohh (3.18) 
IV. Experimental Indications 
ROBERTS [9] concluded from his experiments that, when 
Vik ty. vee Vee, Ae = CONSE, (4.1) 
the relations 
too = bg, hy lee (4.2) 
hold. From (3.6), (3.9), (4.1) and (4.2), 
B2=0, Bs+0. (4.3) 


Hence (3.14) should be satisfied with k=0. Assuming f, is such that (3.11) 
admits suitable solutions, we see from the discussion of Section III that the fluids 
which he observed should be capable of undergoing rectilinear motion through 
tubes of arbitrary shape, to which they adhere. 

From data obtained from different types of measurements, Marxovitz [1] 
concluded that (4.2), hence (4.3), does not hold. The quantities »,, v, and y which 
he uses are related to the quantities B,, 8, and X used here by 


3% = (2B. + Bo) X 
3 y= 2(B, — Be) X, (4.4) 
tye Cs 
As he mentions, data obtained from torsional shearing of fluids between parallel 


circular plates indicates that, for a wide range of values of y, his equation (26) 
holds quite accurately. It is equivalent to 


BX +2 bX (4.5) 


where « is a constant. As is indicated by the data summarized in MARKOVIT2’ 
Figure 5, the normal stresses measured in torsional shearing and coaxial cylinder 
instruments are equal, within experimental error, which means that the left side 
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of MarkKovitz’ equation (18) should be equal to —3(r,+27,). Translating this 
equation into our notation, we obtain 


d 


3 (Bs — 2B.) VX = Bs VX + 2 ana B.X. (4.6) 
From (4.5) and (4.6), 
CS eee ‘ 
oh 
7 20 2B 
B,VX == +, 


where # is a constant. From a more detailed and somewhat different analysis 
of existing data, MARKoviTz concludes that »,=—yg, which is in agreement 
with (4.7) if and only if B=0. This suggests that 


Vy yr 2 
BX Ee , Bs /x oa a7 (4.8) 
in contradiction with (4.3). If (3.14) held, we should have, for some constant , 


By \x= ae (4.9) 


Comparison of (4.8) with experimental data indicates that, if we choose units 
such that /x =y,, where y, is MARKovitTz’ “reduced rate of shear’, we obtain 
good agreement in the range 10°<|/X S104. Within this range, the rather 
meagre, scattered data indicates that 


In B, VX |yx=10 —Inp, /X |yx=10 <1 
If (3.14) held, this statement should hold with f, replaced by f,. According to 
the viscosity data used by MARKOVITZ, 


dinp,VX __ ae 
imyx 7 1,43 —0.12In|/X, 


from which we find that 
In py )/X |yx=10 —In B /X |yx=s0 = 1,01 


approximately. It would thus seem that (3.14) does not hold, hence that the 
fluids observed are incapable of undergoing rectilinear motion through tubes of 
arbitrary shape. 

From (3.5) we see that neither (4.3) nor (4.8) is consistent with the Reiner- 
Rivlin theory. Experimentally, neither (4.3) nor (4.8) has been established 
conclusively and it is uncertain why contradictory conclusions derive from 
different experiments. At least until this discrepancy is resolved, it would be 
premature to conclude the attractively simple Reiner-Rivlin theory to be inade- 
quate. The theory considered in this paper is not as general as that proposed by 
GREEN & RIvLIn [14] and is, in some respects, less general than that of Nort [15], 


The experimental data referred to above is mainly for polyisobutylene solu- 
tions and may not be representative for other types of fluids. 
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V. Flow Through Elliptical Tubes 

In view of the impossibility of conducting experiments using all possible 
shapes of tubes, it is desirable to have some assurance that, if secondary flows 
are not observed in tubes of some particular shape, they will not be observed 
in tubes of any other shape. It seems probable that (3.14) must be satisfied for 
rectilinear motion to be possible in an elliptical tube. This can be proved for the 
class of fluids for which 6, ==const. For these, there is a unique solution of 
(3.11) which vanishes on the ellipse 


x3 + p45 = 9, 
where p and q are non-zero constants with, p?=+-1. It is given by 
1 


far(gtPe—F), 1r=z57(+?%). 


For tity X=} 4 =A 


xX i fi= PARP hes = 1077 (xe “pe x5) 
== Are Pt) Am ar ha ee 
La orl + 2: 
For (3.12) to be satisfied, we must have 


a = 2(Bof,;),;= 2X f 4BldX + 2Bof i; 


= 8r[(1 + p*) X — arp (f+ 19")] dB2/dX + 4Byr(1 + p?). 
Assuming that 7+-0, f and X are functionally independent. We thus have 


(5.1) 


ae = — 327° p2dp,/dX, 
from which 
Ges — 32 pOt+t, a= dX +7, 


where 6, ¢ and 7 are constants. Making these substitutions in (5.1), we obtain 
t= 127(1 + p2) OX — 327° p2q?2d + 4nr(1+ p%), 


which is impossible unless 6=0, in which case B= =const. Thus (3.14) must 
hold for rectilinear motion to be possible. Solutions illustrating secondary flows 
in elliptical tubes are given in [3] and [7]. 

To our knowledge, there is no very conclusive evidence that non-Newtonian 
fluids adhere to solid boundaries. For fluids such that (3.14) does not hold, it 
would be reasonable to expect secondary flows to be observed even if some 
slippage occurs. The effect of slippage should perhaps be considered in deducing 
the functional forms of the 6’s from experimental data, this being one possible 
source of the discrepancies mentioned earlier. 

It is not completely certain that secondary flows will not be observed in the 
steady flow of fluids through tubes when (3.14) is satisfied. The possibility exists 
that rectilinear motion and some more complicated type of flow are both possible 
in some fluids and that these fluids may prefer the more complicated pattern. 
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Asymptotic Solutions 
of ui +B eu" +au' + Bu) =0 for Large |4| 
ALBERT L. RABENSTEIN 


Communicated by C. C. LIN 


1. Introduction 
The asymptotic theory for the Orr-Sommerfeld equation of hydrodynamic 
stability has a large literature. The structure of the equation is such that it 
possesses a turning point; and herein lies the main difficulty in the development 
of such a theory. The equation is of the form 


2 2 
(1.1) G+ BLE, (y) pe +B Only 2) pe” = 0 
n=0 n= 


where the complex parameter A is supposed to be large in magnitude, the func- 
tions P,(y) are analytic with 


(1.1 a) Fy(Yo) =0, P' (Yo) = 0 
and 
(1.1) Only, a) ~ >i ar) | 


where the functions q,,(y) are analytic. In the hydrodynamical case, 7 (y) = 0. 
We are interested in the behavior of solutions of (1.1) in a neighborhood of the 
turning point y) when |A| is large. 

In order to use the method of a comparison equation to obtain such an 
asymptotic theory, one needs an equation which is sufficiently similar to (1.1) 
but which can still be solved explicitly. Such an equation is 


(1.2) wi?t A2(zu"+au'+ Bu) =0. 


Here z is the complex independent variable and « and 8 are complex parameters, 
dependent on / in such a way that they are uniformly bounded for |A|=| A | 
for some Ay. This equation has a turning point at z=0, and is indeed a special 
case of (1.1). Furthermore it can be solved using the method of Laplace integrals, 
since the independent variable appears in only one of the coefficients and appears 
there linearly. 


The complete theory for (1.1), based on the functions defined by (1.2), will 
be presented later*. 


* Preliminary forms of this theory have been reported by Lry. See the Bulletin 
of the Amer. Math. Soc. 62, 394 (1956), and the Proceedings of the IX® Congrés Inter- 
national de Mécanique Appliquée, Université de Bruxelles 1, 137 (1957). 
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The present paper concerns the behavior of the functions defined by (1.2), 
for large |A|. This equation is clearly a generalization of that considered by 
Wasow [1]. It turns out that it is necessary to consider (1.2) with general « 
to make the general theory for (1.1) feasible, even in the case where P,(y) =0. 
Also, there are certain general relations that hold for the functions u(z; A, «, B), 
and these are not evident for the functions u(z; 4, 0, 8) considered by Wasow. 
Indeed, one of the important differences between our theory for (1.1) and those 
of Wasow [2] and LaNncER [3] is the use of the comparison equation (1.2). The 
present theory follows more closely the earlier work of LANGER for the second order 
equation [4]; it is simpler in structure and more readily applicable. 

Series which satisfy (1.2) formally can be found by well known methods. 
These have the forms 


2 91nl2) A ek onl) A ae, 


Siazk * 2iAz8 n 
e Duel) Ls os doa At 


In these, oj,,(z) is an entire function, og9(z) has a branch point at z=0, and 
5 


a 

639 (2) =O49(z) =22 4. Since all solutions of (1.2) are entire functions of z, the 
multi-valued nature of the expressions (1.3) precludes them from representing 
actual solutions of (1.2) in a full neighborhood of z=0. It is known [4] that one 
of these formal series may represent one solution in a certain region of the z 
plane and represent another solution in another region. In any case, the represen- 
tation is known to be valid only for |z| =| z|>0, where 2 is a fixed number, 
so the immediate vicinity of the turning point is not included in the region of 
representation. 


(1.3) 


If we put 7 =A} z, w(z)=7(y) in (1.2) we obtain the equation 
(1.4) ering’ +ay +paty=0. 


Some of (1. 4) having the form of power series in A~# will be valid for |7| <| 7o| 
or |z| <| A|~*| | for an arbitrary but finite y). Thus, for any solution of (1.4), 
the region of validity in the z plane becomes vanishingly small as | A| becomes 
large, while the representation of a solution by one of the expressions (1.3) still 
holds only for |z|=|2 |. Also there remains the problem of connecting the 
solutions of (1.4) with the solutions represented by the expressions (1.3). 

In this article we obtain solutions of (1.2) by Laplace integrals. Certain 
relations which hold among the solutions defined by these integrals are deduced 
with the aid of CaucHy’s theorem. These relations enable one to describe fully 
the asymptotic behavior of each particular solution in a full neighborhood of 
z=0. 

2. General formulae 


Seeking solutions of (1.2) in the form 
u=fe'f()dt 
E 
where ¢ is a complex variable, we find that f satisfies the first order equation 


CS i ay aU 
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where / and C are to be such that [e*’#/]-=0*. Then a solution of (1.2) is 
(2.1) u = ft**exp (4478 —Btt+ 22) dt, 
C 


where C is to be such that 

(2.2) [i exp ($4728 — Bi++ 2t)|c=0. 

If repeated differentiation under the integral sign is permissible, there follows 
from (2.1) that 

(23) ale) yea tn). 


It will be seen that this is the case for every solution considered here. Thus, if 
the asymptotic behavior of a solution (2.1) can be obtained for general «, the 
behavior of its derivatives follows from (2.3). From (1.2) and (2.3) we also obtain 
the recurrence formula 


(2.4) u(z,a+ 4) = — Afzu(z,a + 2)+au(z,0+1) +P u(z, a)]. 


In what follows, we suppose A to be large in magnitude and to have constant 
argument. Having chosen a particular value for arg A, by 4? we shall mean 


A? = exp {fp (log |A| + i arg A)} 


for all complex ~. We assume B=+0, but « and # are otherwise arbitrary. As 
remarked in the introduction, they may depend on 4 in such a way that they 
are uniformly bounded for |A|2|A,|. We take for arg B the value which lies 
in the range 


(2.5) 2argASargB<2a-+ Zargd. 
We now define particular solutions of (1.2) by making appropriate choices 


for the path C of integration in (2.1). Because of the essential singularity of 
the integrand at t=0, paths tending to the origin must do so with direction 


(arg + y=] (mod 2), where —1<y<1. Because of the factor f, the inte- 


grand in (2.1) is in general multiple-valued. The quantity 4? has a negative 
real part in the sectors 


: a 1 2 
T, 6 Nt MR Corts Sy arg A 
/ : S ele 2 786) 2 
2 6 3 ats A<argt< 6 +3 arg 4 
: as = 2 ss 2 
T; 7 bg ate A <argi< 6 1g arg 4. 


These are the unshaded sectors in Figure 1. We also define 


15 — G+ Garg a<argi< — 4% 4S ared. 


* For a discussion of the use of the Laplace integral in solving differential equations, 
see CODDINGTON & LEVINSON, Theory of Ordinary Differential Equations, New York, 
1955, p. 170. 


— ae ee 
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In Figure 1, 7, and 7, have the same location, but the function ¢ will have in 
general* different values in 7, and 73. We define three solutions u=A,, 
k=1, 2,3, by taking paths C(A,) going from infinity to infinity in two unshaded 
sectors. The path C(A,) goes to infinity in J, and J, the path C(A,) goes to 
infinity in 7; and 7,, and the path C(A;) goes to infinity in 7, and 7h. 


Fig. 1. Typical configurations in the ¢-plane 


We obtain three more solutions u= B,, k=1, 2,3, by taking paths from the 
origin to infinity. The paths C(B,) are described as follows. The path C(B,) 
goes to infinity in the sector 7, and to zero with direction arg?—arg f. The path 
C(B,) goes to infinity in J, and to zero with direction argt=argf. The path 
C(B3) goes to infinity in 73 and goes to zero with direction arg¢=arg B — 22. 

It will be convenient to define a solution B,=e~*"'*B, by taking a path 
C(B,) going to infinity in the sector T, and tending to the origin with direction 
arg t=arg B— 22. 

We also define a solution u = B, corresponding to a path C (By) which emerges 
from the origin with direction arg!—arg fB—2z, circles the origin, and tends 
to zero again with direction arg ¢=arg f. 


* In other words, J and Z, will lie on different sheets of the Riemann surface 
for the function ¢* = exp [« (log |¢| +7 arg ¢)]. 
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Convergence at the origin is not altered if the direction of approach is changed 
by an angle less than }z, in accordance with the remarks following (2.5). Use 
will be made of this later. For each path C described above, the integral (2.1) 
converges uniformly for z in any finite neighborhood S of z=0 and for « in any 
bounded region of the « plane. Also, the condition (2.2) is satisfied. Then dif- 
ferentiation under the integral sign is valid, and (2.3) follows. These paths of 
integration are shown in Figure 1 for various values of arg f. In the visible 
portion of the Reimann surface, argt ranges from arg f—2z to argf. The 
dotted parts of the paths are the parts which lie on other sheets of the surface. 
The paths which go to infinity do so with directions independent of f. At the 
origin, the direction of approach varies continuously with arg f, for arg f in the 
range (2.5). As a consequence of Caucuy’s theorem, the following relations 
hold among the solutions: 


B,+ 4,= Bs, 
(2.6) Bec As =D, Ba (B, = B, e~?7'4), 
B,+ A,= Bs. 


The paths of integration have been chosen so that these relations hold for all 
values of B. 

We now seek the asymptotic description of each of the functions defined 
above in a full neighborhood of the turning point z =0. The development follows 
closely that used by Wasow in [J]. 


3. The solution B, 
We have 


(G1) By= f @#-*exp (zt — Bi>) exp (22-28) dt. 
C (Bo) 


poly)  a 


converges uniformly for ¢ in any bounded domain and |A|=|A,|>0, and since 
the rest of the integrand in (3.1) is uniformly bounded for ¢ on C(B,), termwise 
integration is permissible and we obtain 


Since the series 


(3.2) Bo = 09(2, a) + 3A? bo (z,0% + 3) +--- 

where 

(3.3) boo) = | PA expat — Biayas. 
C (Bo) 


This may be expressed in terms of a Bessel function. The change of variable 
t=e ™' Br gives 
Dole; ce) eter as Bas tes *exp (z - Pe) ar. 
Tt 
C’ (By) 
Comparison with the integral representation formulas for a Bessel function of 
the first kind ([6], p. 176) gives 


(3.4) bo (2, %) = — 2rre **% (2b B-H1-o (2 28 Bh). 
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This may be written 


(3.5) bg (2, «) = — 2aie—™** fr oy eS : 


Thus By tends to a single-valued analytic function of z in a full neighborhood 
of z=0. We note here that 6)(z,«) satisfies the second order equation 


(3.6) zu’ +au'+ Bu=o0 
obtained from (2.1) by formally allowing 4 to become infinite. 


4. The solutions A,, 
We have 


(4.1) Ap [exp (E2728 Set — 8 i at 
C (Ax) 


We shall evaluate A, asymptotically by the method of steepest descents ([7], 
p. 503). For this purpose let us write 


(4.2) Hr 
and 
gA?7B 4+ ct =FAPVG+tyzr=u(se—r). 
Then 
(4.3) B= AZ = — ye, 
whence 
(4.4) y= — At, w=tA(2), 


where the same branch of the square root is to be taken. From the formal solu- 
tions (1.3) we expect the quantity 


(4.5) E=2i,A=24u 
to be important. In terms of €, (4.3) gives 
(4.6) v= MBE, z=—A FSP )?, 


where the same branch of the cube root is taken. The proper branches will be 
determined in the sequel. Let us now write 


pA 28 + 2t=E(Gt° — 37) = |é|F(z) 


where 

(4.7) faj=er(Gr Fr), od =argé. 
Then (4.1) becomes 

(4.8) Ap =A ( eae eG) Aye 2 
where 


F(G8 2) =e texp{— BANE) A, 


f(t) is given by (4.7) and C’(A,) is the image of C(A,) in the t plane. Clearly 
the location of C’(A,) depends on the choice of a branch for y in (4.6). 
Arch. Rational Mech. Anal., Vol. 1 29 
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We consider first the case where the variable & is large in magnitude, say | 
|E|=|&| for some &. | 
The function Re{f(z)} has a saddle point at t=1 where /’(t) =O, and 


(4.9) (Sd Seep ae 
The path of steepest descent through t=1 is given by 
(4.10) Im {f(z) — f(1)} = 0. 


Fig. 2. Path of steepest descent 


t-plane 


m . 1 ; » - 4 
Y=>5 0] 1) Se ea j= 1, 2,3, 5,6, 7—2=<— 0—-2m2n= nH, y=argt 


To plot this path, put t—1=ge'®. Then (4.10) becomes 
o?[o sin (3@ + 9) + 3 sin (20 + #)| =0. 


For # in each of the ranges 
(4.44) —na+t+2mn<0<at2mn (je 0, tale Los eG 


the path of steepest descent is as shown in Figure 2. In the Figure the unshaded 
sectors are those in which e'” 73 has a negative real part. For each range of @ in 
(4.11), the path of steepest descent goes to co in the two unshaded sectors indica- 
ted. As # passes out of any one of the ranges (4.11), the path changes abruptly 
and goes to oo in the third unshaded sector. 

The idea now is to choose, for each A,, a branch for y in (4.6) so that for & 
in the region specified by (4.11), C’(A;) goes to infinity in the same two sectors 
of the t-plane as does the path of steepest descent. Then we can take for C’(A,) 


this path of steepest descent. This requirement is satisfied if we take m =k —2 
in (4.11) so that 


(4.12) — a+ 2u(k —2)<0<a+ 2n(k — 2) 
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ut 


and take 
(4.13) y = AB|Z El eh? 


for 9 in the range (4.12). To see that this is so, let C(A,) go to infinity in the 
t-plane with the directions 


argi= 42 ss arg A + Sl (k — 2). 
3 3 3 
Then since arg t=arg ¢—argy, we see that C’(A,) goes to infinity in the t- 
plane with directions 
argt= + if aa Bic 20 ae) 
re) e) 


Inspection of Figure 2 with m replaced by k—2 shows that the choices (4.12) 
and (4.13) are satisfactory. 


To evaluate A; by the method of steepest descents, set 


(4.14) Hy fd Senta ee r= 18 a se 

Choose 

(4.15) s = (— 1)** 17 e8*9 (¢ — 1) (r + 2)3 

with the branch cut along the negative real axis from t = — 2 to t= ov, and take 


the branch of the square root which is real and positive for t real and t> —2. 
Then as t traverses the path of steepest descent in the direction corresponding 
to that for C’(A,;), s increases through real values. The derivative ds/dr is finite 
and non-zero in the cut plane, so near s =O, Tt is an analytic function of s. For 
real s, and |s| = 6>0, 


F(x(s),€, A)| S K,| st-?)|, 


for real s. If we expand F dt/ds in a Taylor series about s=0, it follows from 
Watson’s lemma ({7], p.501) that the series obtained formally by termwise 
integration over (— oo, co) represents A, asymptotically for 


(4.16) — 3% + 2n(k—2)<0<2% + 2m (b— 2). 


Watson’s lemma extends the range of # in (4.12) by $z in each direction. It 
can be verified from (4.4) that taking # in the range (4.16) and choosing y as in 
(4.13) is equivalent to restricting y=arg< to the range 


(4.17) Bp hy are ay on a Gee + args 
3 2) 

and taking ; 

(4.18) ye * A|z|Per? 

(4.19) E=2iAlz|P et? 


29* 
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The results of the integration are 


(4.20) Api Yr(— 1)* AOD EVIOH e#[1 + (EB 2 — le — 3) + EF $d, 


where #=arg é is in the range (4.16). This can be written 


3 
2 


Satie fk: 


(4.21) A,=iVa(—1)ke 2 3° 


x ze 2 ide 
zg? 4 e 2 


3 2) 4 Th 403] 


where y = arg z is in the range (4.17). 

Let us consider again the function é in (4.5). In the z plane, Re(€) is zero on 
three rays, denoted by C,, k=1, 2,3, extending from the origin. These rays 
are given (mod 22) by 


2 2 2 
Ci: are = Se: : ee 3 arg A. 


These divide the finite neighborhood S of z=0 into three sectors S, of angle $2. 
We take S, to be the sector opposite C, and take S; to be closed. It can now be 
verified that the choice of a branch for & in (4.17), 
(4.19) is equivalent to restricting z to S—C, and 
choosing the branch of € for which Re(é) 20 in 
S, and Re(é) <0 in the other two sectors. 

If we examine the relations (4.20) we see that 
for zin any closed subset of the interior of S,, A, 
tends exponentially to zero as A becomes infinite 
and becomes infinite exponentially for z in a 
closed subset of the interior of either of the other 
en two sectors. The behavior is described by saying 
that A, is subdominant or recessive in S, and 
dominant in the other two sectors. 


C\ 


Cc 


Nn 


Fig. 3. Diagram for choice of branches 


The behavior of A, near C, can be obtained from the relations (2.6), once 
the asymptotic properties of the solutions B, have been investigated. 

For the case when |é|<|&| we put ¢=A5o, z=A7'y in (4.1). We note that 
|&|<|& | implies |7|<|mo|. These substitutions give 


(4.22) A, == At tae i ope exp(& +7 o| exp(— BA-to) do. 
C*(Ax) 
If we write 


n=0 


oO 


and integrate term by term, this being valid here, we get 


en co (4 yeh 
(4.23) Aya BOD YE ala — 2), 
where 
(4.24) x(n, &) = [ ox *exp(2 +0) do 
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We observe from (4.24) that 


qd” 
(4.25) Cre) — ey, a+n). 


The functions g, (7, «) satisfy the equation 
(4.26) x? tne’ +ay’'=0 
obtained from (1.4) by putting formally A= oo. Then because of (4.25), 


(4.27) Ban, % + 4) +7 &4(y,% + 2) + a8, (9,0 + 1) = 0. 
This holds for all «, and therefore 

(4.28) 8a ly. % + 3) + 79 8a (y,% +1) + (% — 1) ge (y, 4) = 0. 
Hence g,(n,«) actually satisfies the third order equation 

(4.29) KX +nx +@—1)%=0. 


For a discussion of this equation see LANGER [8] and HERSHENOV*. 


5. The solutions B;,, k = 1, 2, 3 


We have 
(5.4) = hfettrexp HAP Bee 4— Bi) ai 
6 (Bx) 
eet 
(5.2) by (zen) = fc *expltz — ft) dt 
C (Bi) 


This integral converges for all « if Re(tz)<0 for large |¢|. The requirements 
Re (tz) <0 and Re(A~*#5) <0 for large |¢| can be satisfied by deformation of the 
path C(B,) if and only if z€ S—S,. Now 


pi atte S Le ne 78 an JA-2+) {3 ( (n+1) hy (t, A), 
where h,, (¢, A) is analytic in ¢ for fixed 4 and is uniformly bounded for ¢ on C(B,) 
and |A|=|A,|>0. Therefore (5.1) gives 


(5-3) By = 6, (2, %) + --- aT (4 4-2) b,(z,a + 3n)+ R,, 
where 
(5.4) RHA Pern f ptsery—2 exp (tz — Bi) h,(¢, A) dt 
C (Bx) 
Te }2) | 2,| > 0, then 
(5.5) OCA al): 


In the case where | &|=|&| but z is not bounded away from the origin, put t=¢2 
n (5.4). Then 


fi, 26 4) f= Ga sn > SER es MLE \ip® 
(5.6) Vea | z ike exp(r = )m (eA) ee 


C’(Bx) 


* Master’s thesis, Department of Mathematics, M.I.T. (1957). 
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For n> —2—1Re(a), the integral here is uniformly bounded and 


(5.7) R,, = Oilgh= are igri 
Later it will be seen that regardless of the value of « 
(5.8) B, = 5, (z,«) [4 + O0(E°)]. 


For some values of « the error term O(&*) can be replaced by one of still smaller 
order of magnitude. 

To express the functions 6,(z,«) in terms of special functions, we make the 
change of variable 


(5.9) s= tz prt 
n (5.2), obtaining 
= (pry [otexp 


C*(Br) 


zt pi(s — >) ds. 


S| 


Comparison with the integral representation formulae for the Hankel functions 
(see [6], p.179) shows that with the proper choice of branches in (5.9), 


(5.104) by = ni (Bt ah) HY, (0p! 2), 
(5.10b) b, = mi (B-* z4)1—* HY, (2B 28), 
(5.10c) bs = 21 (e?** B—-# 2h) -* A (2,88 24) 


where 2? =| z|? e*?, 


(5.11) he oily (3 2) ae argd<g< 7% EaN (k — 2) ? arg, 
3 3 3 3 3 3 


and arg B+! = +4 arg B, with arg B as in (2.5). Examination of these functions 
for various values of « proves (5.8). We remark that the functions 6,(z, «) satisfy 
the second order equation (3.6). 


We now give a more detailed analysis for the case k= 2, using the formula 
In) pai79. For a0 |t] we wish to have 


2B + arg a <argi< + Sarg d 
according to the choice for C(B,), and we want 
> + 2n 2a < arg (tz) << i +22 


for some integer 1, so that the integral b, converges. These will hold if we take 


argi = —qargz+ a+ argi+ ina, 
where 


Fmt ann — = arg d< args < 2% + 2nx 5 arg A, 


for large |¢|. For the moment take arg 6! = 3 arg B + ma. In the s-plane, the path 
of integration tends to infinity in the ar ceien 


args = a+ Zargz+ gargd+inn—ma—targp=n—ow. 
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(See [6], p. 179 for the significance of w.) The direction ue approach for C(B,) 


at the origin in the ¢-plane can be taken as arg t=argf +95 , where —1<y<1. 
We choose y so that at the origin in the s-plane 


args = | argz | + ar p-ma+y—=o 
3 Sasa al bees ‘ 
The two expressions for w give 


VIC 


zy = 1 3 
y= = (2m Z arg A i arg 2). 


To have |y|<1 we must have m—n=0. Next, we find 


4 14 
ig oN — arg (6) 4) << — 2a, 


so we take n»=0, and then m=O. Here =F aO= La. The formula in [6] is 


given for |w|<z, but extending the range of w simply gives an analytic continua- 
tion of the Hankel function. The relations (5.10b) and (5.11) with k=2 now 
follow. 


The asymptotic behavior of B, in S, can be obtained from the relations 
(2.6). We note that B, is dominant in S,. For z in S—S, and |z|2|z|>0, 
B, tends to a finite non-zero limit as A becomes infinite. We say that B, is 
balanced in S — Sy. 

We now examine the functions B, in the case where |&|<]é|. Here it is 
impossible to proceed as with the solutions A, because the path of integration 
goes to the origin. 

When « is not an integer we can solve the equations (2.6) for the three functions 
B, in terms of the functions By and A,. We get 


—Onias 
Ue AS een neers 


(5.124) B= —Quia , 
Le 
(5.12b) B= A, +4, vars ‘ 
1 — eet 
(5.42¢) peu iae aii neem 


—2n1 
Naar m1% 


The asymptotic expansions for the functions By and A, when | &|<|&| are given 
by (3.2) and (4.23). 


Since all solutions here are entire functions of «, we can pass to the limit 
as « tends to an integer m, obtaining 


(5.13) — ee (A, 4 Aa +A,— By) — 271 A,| SS 
[ a 

Ouge) Ba= sri fgy Mat dat 4s Bo], 

(5.13) Bees 3a — (A, 1 A, Tale Bo) 2mi Ay) 
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It is now necessary to examine the quantity 


(5.14) W,= [5 (Ar + 42+ 4s— Bo), 


Let K =C(A,) + C(A,) +C(4,), P=K —C(B,). Then 


6 —2 { 9-248 = 
W, =\— |? “*exp|— 77°F -2t— pt i ? 


A fae log texp(= AB 4+ zt pe) dt. 
tig 
If we put 
(5.15) b= z= i", 
this becomes 


W,, = Ae 5 (log A) f o7texp(= Lo B )do | 


3 Bo 


[T* 


a o”* log o exp(5- + 0 — 124]: 
* 3 o 


P 
The first integral here vanishes because of CAUCHY’s theorem. Then 
[oe 


W,, = 8O-D | as a (Fy ii o”* logo exp (7 +n o| do — 
BE ! 5 Pt 


aos 
J 


* ie o"* Jog o exp (= +Ho—- E Jao : 
3 136 
C*(Bo) 
Reversing the transformation (5.15) in the last integral here, we find its con- 
tribution to be 


S (log A) ! i* exp(5 A? B+ zt Br) at— 


C( 


— f #* log texp(4 A2B+ zt pr) dt 


5 
C (Bo) 


te) 
i en (n, % — m|, = 


(5.16) eR 


The functions g,(7,«) are defined by (4.24). The asymptotic expansions of 
the functions 6, when « is an integer m are then given by (5.16) and (5.13) 


By taking «=0, B=1, these results may be shown to agree with those of 
Wasow (([2], p. 227, equations (2.8a,b,c,d)). The relations connecting his 
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solutions U,, U,, U, with those considered here are 


(pel, 
Oe BE 
OS ore 


6. A set of independent solutions 
We consider here the set of solutions B,, B,, A,, Ay. Similar results will hold 
for any set By, B;, A;, A,, where no two of 7, 7, k are equal. The solution By is 
balanced everywhere in S. The solution B, is balanced in S — S, and is dominant 
in S;. The solution A; is subdominant in S; and is dominant in the other two 


sectors. This set will be shown to be linearly independent for || sufficiently 
large. 


a C(a) = tems (ao 
Then for z in S— S, and |£|=|£ | we have 
(6.1) By = bo(z, a) + O(a), 
(6.2) Ba= 0, (eee) dt 0 (S=*)/\; 
(6.3) ee ee eo. 
(6.4) Aya —ie* Cla) A278 22 8 Tt LOE, 
where in (6.2), (6.3), (6.4) 

4m p 


3 


Let W(f,,...,7,) be the Wronskian of 7,,..., 7, with respect to z. Using (2.3) 
to evaluate the derivatives of the above functions, by calculation we find 


W(Bo, By, Ay, Ag) = 271 29"? 2 W(bo, bs) [1 + OE], 
An? Care Be A eto Oe 3). 


— =* Sarg <argz < —Zarga. 
3) 


(6.5) 


In the case where | &|<|& | we have 


(6.6) B= BOE ets a a OM A= 3), 

(6.7) A; = AY [95 (y, «) + O(A-4)] Get) 
Using (5.12c) we find 

(6.8) W(Bo, Bs, Ay, A») ae eSEr W(Boy, 41, Ae, As) ; 


In the case where « is an integer 7, we must take the limiting value of this ex- 
pression as « tends to ». By use of (6.6) and (6.7), this becomes 


* px—1 —382140 
— 221 p* e 


(6.9) I'(a) 1 gente 


a+? W (gt, 2, 83) [1+ 0(A-4)], 
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where W, denotes the Wronskian with respect to the variable 7 and the primes 
denote differentiation with respect to 7. With the aid of (4.29) there follows 


(6.10) W, (1, 82, 3) = (1 — «) W, (61, 82, 83) - 


Since the Wronskian on the right is independent of 7, it may be evaluated at 
n=0. The quantities g,(0,«) may be expressed in terms of gamma functions 
((9}r p,225).. Ineiaet, 


oF 3 am (a—1) (k—2) 
é 
g, (0, a) a fw ( Ona = phy A= pig essa ° 
C*(Ax) bey 


By the identity ([9], p. 225) 


ieee Pe = at = e) Avs ee es ot) 


and (4.25) there follows 


set 2 

(6.11) WA (61> 82» 83) = Tee . 

From (6.9), (6.10), (6.11) and the familiar identity 
I'(2—«) (a) = ANE?) 


sin 7 a 
we obtain finally 


(6.12) WB By Ay ApH ane t Pa ew Oe 
Then from (6.5) and (6.12) we see that if |A| is sufficiently large, the solutions 
By, Bs, Ay, Ag are linearly independent. 


A summary of the asymptotic properties of these functions is given below. 
The behavior of the other functions may be calculated from this information by 
using the connection formulas (2.6). 


(i) The case |é|=|&|. 


Let €= 3742, and in the formulae below let arg z lie in the appropriate part 
of the range 


2 2% 
— 2m — = arg A< arg z<—- — 7 arg A. 
2) J 


(6.13) Bo= bo (2, a) 2:0)" for zan S, 
h 
ba bo (z, a) = — 2uie*i(B-4 atyd—2 | (288 2A), 
PES Tone: 
(6.14) A,=—C(a)A 2 22 4e* [1+ 0(€>)] for z in S—C,, 


3 (of 


(6.15) Ap=—td 2 e**C(q) z? tag [1+ 0(€>)] for sane SG 
(6.16) By=mi(e?** Bo) 2h)!-* HO (288 A) [4 + O(E*)] for z in S—Sy, 
(6.17) Bs=—C(a)A. 22! 7 e*(1+0(€4)] + 

+ wie Ona (Bodine HY (264 A) [4 LO(E?)] for z in S,—C,. 
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The expressions for A, and B; near C, and for A, near C, are more complicated 
but may be found easily from the relations 

B,= By + B,— Ay, 
(6.18) Ay Bei See by Ay 6 A 

Ay Bye) Bae A, Ay. 
These follow from (2.6). 

(it) The case |E|<|&p| (for z in S). 

(O19) BB y= bye) 420.4) S207 2 '* pe 
(6.20) A;= MO” [e;(y, «) + O(A-4)], 


1 Bie 7k (tte 
ag ee 


—nmi4 


(6.21) By = =F {AH (#4 + ge + g9) + O(A-8)] — do(z,«) + O(A)} 


when « is not an integer. 


2\(n— . a) , | —% { 
B, = [eaten + > ae | + O(A-#)h + 


k=1 


(6.22) 1 ; 1 [ @bo(z, «) : 
+z; (log 4) [bo(e,m) +0(a-y]— 4 [PE | + os) 
when « is an integer 1. 
In these formulae . 
s(ya) = f ot *exp(2 +09) do, 
C*(Ax) | 
fe re (4—1) (k—2) 
8x(0,%) = 21 gy 
is 
For z=0, these formulae (27) yield 
- IK Ra— 1 
(6.23) By(0) ~— 2a7e re F(a)’ 
ae ge Ne 
(6.24) A, (0) ~ 221 2 eae Ae AY, 
ere 
2) 
= rs aE Coss = tae 2) 
2 Sa 3 ae 
ee eae ae AB (a — 1) + 
(6.25) a fame 
: T' (a) sinza eam) 
n—4 
of 3 
Bz (0) ~ Ab Fira) x 
3 
CATO ery ee, 1 1 4—n 
(6.26) x |anie? a eae eines (it 1) + 3 log3 + 3 | 3 ] 
2 at —1 1 a —ni4 CP Noe al == 
je Br Fay 084+ ae aerate (a = n) 
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(11) Bounds. 


Let &; be the branch of = 27,25 which has a positive real part in S; and a 
negative real part in the other two sectors. Let E = E(z, A) be a generic symbol 
for a uniformly bounded function of z and 4. Then the following relations hold. 


For |&|2|& | and z€S, 

(6.27 a) Bota 
For |&|<|&| and z€S, 

(6.27b) el 


For |é|2|&5), 


(6.28a) A, 4s Es 6H, 2€S—C,; =1,2,3). 


For-zeCG;, 
1D, Rew 45.0 /5= 0, et 28 
a Seen ee (eee) = og ia 
(6.28b) A,;e°"i =~ A3 ee BE tes Rees oy Or c= 0, — 1, ee 
We EOS ere Rea >3 
For |é|<|&| and z€S, 
(6.28c) Agente Atte. 


Hor |élz\2,| and 26S @ =f) 24 jeer tag 


18, ReGeSA oes 
2 2 
i ee Eee eee 
(G20) Geet ae ea ae 
ear Rea =3 
(log ) F, a= 1. 
For |é|<|& | and z€S, 
E, Rea<1, +1 
(6.29b) B, et ReGe+*d — | (log 2) E, a=1 


APO=DIEY « URGE ose 
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Note added in proof. Professor TRUESDELL has kindly pointed out the connection 
between equation (2.3) and his own work on the equation He (25 0.) == shi (epoca tod) 
iB 


(Annals of Mathematics Studies, Number 18, Princeton University Press, 1948.) Use 


of the theory for the F-equation would have shortened several of the formal steps 
in this paper. 
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Stirungsrechnung 


hei Eigenwert- und Verzweigungsaufgaben 


JOHANN SCHRODER 


Vorgelegt von L. COLLATZ 


Eine gegebene Eigenwert- oder Verzweigungsaufgabe kann man oft naherungs- 
weise lésen, indem man von einem ,,benachbarten‘‘ geschlossen lésbaren Eigen- 
wertproblem ausgeht. Um verschiedene Aufgabentypen (Matrizengleichungen, 
Differentialgleichungen, ...) gleichzeitig behandeln zu kénnen, verwenden wir bei 
der Formulierung des allgemeinen Problems und in den Beweisen funktional- 
analytische Begriffe und Methoden. Fiir bestimmte (in den §§ 8 und 9) als Bei- 
spiele behandelte Stérungsprobleme bei Matrizen- und Differentialgleichungen 
werden die Ergebnisse dann jedoch so zusammengefaBt, daB man sie unmittelbar, 
also ohne funktionalanalytische Formulierung der Aufgabe, anwenden kann. 


Wir behandeln Aufgaben der Form 


Ayp — Ap = Bid), (0.1) 
lo=1 (0.2) 


in einem Banachschen Raum. Gesucht sind also eine Zahl A und ein Element ¢, 
welche diesen Gleichungen geniigen. A, sei ein linearer Operator, / ein beschrank- 
tes lineares Funktional. Die ,,benachbarte“ Aufgabe 4)~9—Agp=0, lp=1 be- 
sitze einen einfachen Eigenwert A=) zum Eigenelement p=yy. Ist B(A, g) 
linear in g, so ist auch (0.1), (0.2) ein Eigenwertproblem, andernfalls liegt eine 
sogenannte Verzweigungsaufgabe vor. 

Bei der itiblichen Stérungsrechnung fiihrt man einen Stérungsparameter ¢ ein, 
indem man 6 durch ¢ B ersetzt, und berechnet 4 und m dann als Potenzreihen in ¢ 


A=AK(8) = Mot eft eust-, gp=ple)=Ytemteyet:. (0.3) 


Diese Methode wird fiir Eigenwertaufgaben z.B. in [2] und [3] beschrieben und 
in [6] auf Eigenwertaufgaben sowie in [7] (in etwas anderer Form) auch auf nicht- 
lineare Probleme angewendet. Wir benutzen hier stattdessen ein Iterationsver- 
fahren. Das hat mancherlei Vorteile fiir die numerische Rechnung [z.B.: ein- 
fachere und bei jedem Schritt gleiche Rechenvorschrift, Selbstkorrektur des Ver- 
fahrens, weitgehend freie Wahl der Ausgangsnaherungen (siehe Beispiel 1), Kon- 
vergenz auch bei nichtkonvergenten Potenzreihen (0.3) (Beispiel 2) ]. 
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Dieser Aufsatz schlieBt an eine kiirzlich erschienene Arbeit iiber das Newton- 
sche Verfahren [16] an, in welcher auch Iterationsverfahren der Form 


Uy Mae [Go Coen G Uy, (1 SUN PAE .) (0.4) 
zur Losung einer Gleichung 
Gu=0 (0.5) 


mit einem G ,,benachbarten“ Operator G, in einem abstrakten Raum t behandelt 
werden. Von diesem Raum wird darin vorausgesetzt, daB jedem seiner Ele- 
mente w als verallgemeinerte Norm ein Element »[w] eines halbgeordneten linearen 
Raumes %t zugeordnet ist. Konvergenzaussagen und Fehlerabschatzungen fiir 
das Verfahren (0.4) werden dabei durch Vergleich mit einem Iterationsverfahren 
entsprechender Form 


nti =n — [Goel]? Ge, (m= 0,1, 2,...) (0.6) 

zur Lésung einer Gleichung 
G 0.0 (0.7) 

im Normenraum t erhalten. 

Nach genauer Formulierung der Aufgabe (in §1) definieren wir hier (in § 2) 
die Raume i und %, schreiben (in § 3) die Aufgabe (0.1), (0.2) in der Form (0.5) 
und leiten aus (0.4) das Iterationsverfahren her. Dann werden (in § 4) gewisse 
Majoranten ermittelt, welche in die anschlieBend (in § 5) untersuchte Vergleichs- 
aufgabe (0.7) und das Vergleichsverfahren (0.6) eingehen. Um die Konvergenz 
des Verfahrens (0.6) zu beweisen, braucht man dabei nur zwei nichtlineare Un- 
gleichungen fiir zwei reelle Unbekannte zu losen. Das Ergebnis wird in § 6 zu- 
sammengefaBt. Besonders einfach wird der Rechengang, wenn A, ein selbst- 
adjungierter Operator eines Hilbertschen Raumes ist. 


Etwas ausfiithrlicher behandeln wir (in § 7) noch den Spezialfall 


(Ag + 41) p =A®q, (0.8) 
Leal (0.9) 


der Aufgabe (0.1), (0.2) in einem Hilbertschen Raum mit selbstadjungiertem 
Operator Ag, linearem Operator A, und dem Funktional /¢~ = (g, pp»). Die Lésun- 
gen der genannten zwei Ungleichungen (bzw. Gleichungen) lassen sich hier (wie 
auch in anderen Fallen) in geschlossener Form angeben. Die Ergebnisse fiir 
diesen Fall enthalten u.a. auch die Fehlerabschatzungen fiir die Naherungen 
Mot /4 und yy fiir A bzw. ¢, welche in [13] auf anderem Wege bewiesen wurden, 
und zwar einschlieBlich einer Verbesserung, auf deren Méglichkeit F. W. SCHAFKE 
[12] hinwies (s. auch Beispiel 1). 

Der Satz aus §6 wird dann auf bestimmte Klassen von Eigenwertaufgaben 
bei Matrizen (in § 8) und gewohnlichen Differentialgleichungen zweiter Ordnung 
(in §9) angewendet. Unter anderem werden fiir Matrizeneigenwertaufgaben 
Formeln angegeben, mit deren Hilfe man den Einflu8 von Abrundungsfehlern 
auf den Eigenwert abschatzen kann. — Die Genauigkeit der Abschatzungen priifen 
wir (in §10) numerisch an zwei einfachen Matrizeneigenwertaufgaben und er- 
lautern am Beispiel eines Knickproblems (in § 11) ausfiihrlich die Anwendung 
des Satzes aus § 6 auf nichtlineare Probleme. 
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Ein wesentliches Ziel dieser Arbeit ist es, fiir die betrachteten Probleme 
numerisch brauchbare Fehlerabschaitzungen zu bekommen. Das Literaturver- 
zeichnis nennt einige Arbeiten aus der Literatur iiber die Stérungsrechnung, 
welche mit diesem Aufsatz in engerem Zusammenhang stehen und insbesondere 
auch Fehlerabschatzungen enthalten. 


F. ReLiicu [8] untersuchte u.a. Stérungsprobleme der Form 


(A,+¢e4,+ @A,+::)p=Alp (0.10) 


in einem Hilbertschen Raum mit (linearen) Hermiteschen Operatoren A; und 
selbstadjungiertem Operator A, und leitete dabei auch Fehlerabschatzungen fiir 
die Potenzreihen (0.3) her. B. v. Sz. Nacy [/7] verbesserte diese Fehlerabschat- 
zungen mit Hilfe einer andersartigen Beweismethode und erhielt entsprechende 
Ergebnisse auch fiir Probleme der Form (0.10) in einem Banachschen Raum [18]. 
Fiir Probleme der speziellen, aber praktisch wichtigsten Gestalt (4) +¢4,)@ 
=Ag in einem Hilbertschen Raum mit selbstadjungiertem Operator Ay, aber 
nicht notwendig Hermiteschem Operator A, wurden in [/3| weiter verbesserte 
Fehlerabschatzungen fiir die Potenzreihen (0.3) angegeben. Von dem Operator A, 
wird darin verlangt, daB eine Ungleichung 


| Ar el S A ( , Ao@l) (0.11) 


mit geeigneter Funktion Q,(s, ¢) gilt. Dies ist eine Verallgemeinerung der bei 
RELLICH und Sz. Nacy benutzten Abschatzung ||A, || < a||q||+0||4o¢|l. 

P. ROSENBLOOM untersuchte in [9] Iterationsverfahren zur Lésung von Auf- 
gaben der Gestalt (0.1), (0.2) in einem Banachschen Raum mit von A unabhangiger 
Funktion B und anschlieBend zusammen mit H. Brock [1] auch solche Auf- 
gaben, bei denen B von 4 und g abhangt!. Gefordert wird dabei u.a., daB B 
auf einer bestimmten Teilmenge Lipschitzbedingungen 


|B (4,9) -—Blu PI|SP\A—p|, BA 9) -B Aw SElle—yll (0.12) 


geniigt und / Eigenelement des adjungierten Operators Aj ist. Das behandelte 
Iterationsverfahren wird als Verfahren der Form 


Uii= Tu, (Gia Oe ee) (0.13) 


gedeutet und darauf der bekannte Fixpunktsatz fiir kontrahierende Abbildungen 
angewendet. 


Unabhangig davon wurden in [15] Iterationsverfahren zur Lésung von Eigen- 
wertaufgaben der Gestalt (0.8), (0.9) in einem Banachschen Raum untersucht, 
und zwar als Beispiel fiir die Anwendung allgemeinerer und genauerer Konver- 
genzkriterien fiir Iterationsverfahren der Art (0.13). Diese Untersuchungen 
wurden in [/6] (siehe insbesondere § 6: Stérungsrechnung) fortgesetzt, und die 
dort erhaltenen Ergebnisse werden hier nun angewendet. Dabei brauchen die 
Bedingungen (0.12), welche z.B. fiir die in § 14 behandelte Aufgabe nicht erfiillt 
sind, nicht gefordert zu werden, und das Funktional / ist weitgehend frei wahlbar 


1 Bock und RosENBLOoM untersuchten auch den Fall, da8 A, nichtlinear ist. 
Sie benutzten andere Bezeichnungen, die von ihnen behandelten Droplets lassen 
sich jedoch in der Form (0.1), (0.2) "schreiben, 
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(siehe ebenfalls § 11). B(A, gy) darf hier die Form einer Funktion C (A, vy, Ay @, 419, 
..., A,g) haben, welche in gewissem Sinne von einer reellwertigen Funktion C 
majorisiert wird. Die A; bedeuten darin lineare Operatoren, fiir welche Unglei- 
chungen der Art (0.11) gelten. Die hergeleiteten Fehlerabschatzungen sind ein- 
facher und genauer als die in [/] angegebenen. Das hangt damit zusammen, 
da8 der in (0.13) verwendete Operator T wegen des Gliedes A in (0.1) auch bei 
linearen Eigenwertaufgaben nichtlinear ist und eine Lipschitzbedingung fiir den 
Operator T immer nur in einem gewissen Teilbereich erfiillt ist. 

F. W. ScHAFKE behandelte in [12] (im AnschluB an friihere Arbeiten [10], [17] 
liber speziellere Problemklassen) die Aufgabe (0.10) in einem Hilbertschen Raum 
mit selbstadjungiertem Operator A, und nicht notwendig Hermiteschen Opera- 
toren A,, A,,..., fiir welche Abschatzungen der Art (0.11) gelten. Er leitete 
darin u.a. eine einfache, in gewissem Sinne nicht zu verbessernde Schranke fiir 
den Konvergenzradius der Reihen (0.3) her und verbesserte die in [13] bewie- 
senen Fehlerabschatzungen fiir w+, und yy (siehe dazu auch § 7). Im Beweis 
benutzte ScHAFKE funktionentheoretische Uberlegungen und den Fixpunktsatz 
fiir kontrahierende Abbildungen (in anderer Weise als BLocK und ROSENBLOOM). 

Weitere Arbeiten tiber Stérungsprobleme der Art (0.10) stammen z.B. von 
T. Kato [5] und F. Worr [19]. 


§1. Aufgabenstellung 


1.1. Es sei S ein reeller oder komplexer Banachscher Raum, dessen Elemente 
mit @, p, ... bezeichnet werden. ||q|| bedeute die Norm von @. Ay sei ein linearer 
Operator, welcher eine lineare Mannigfaltigkeit Uc G in © abbildet. Er besitze 
den Eigenwert ww) zum Eigenelement wo: 


Ag Vo = Io Yo- 


Die weiteren Voraussetzungen fiir Ay werden zunachst fiir den wichtigsten 
Spezialfall und erst anschlieBend fiir den allgemeinen Fall formuliert. 


I. Spezialfall?. &S sei ein Hilbertscher Raum mit dem inneren Produkt (¢, y) 
und der Norm ||||=V(¢, ¢). Ao sei selbstadjungiert, ji) ein einfacher Eigenwert 
von A, und yw, durch 

(Yo Yo) = 1 


normiert. Das Spektrum von Ay, sei leer in einem Intervall 
Mo —€<A<potd mit: 0<d<oo 


abgesehen von A= [y. 


II. Allgemeiner Fall. wy sei so bestimmt, daB 
Ly = 1 (1.1) 


gilt mit einem (der gegebenen Aufgabe entsprechend zu wahlenden) auf © er- 
klarten beschrankten linearen Funktional /. 


2 Dieser Fall wird im folgenden als ,,Spezialfall 1“ zitiert. 
Arch. Rational Mech. Anal., Vol. 1 30 
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Es gebe ein auf © definiertes beschranktes lineares Funktional ® derart, 
daB die Gleichung 
(Ay— Mol) p=y¥ (1.2) 


(I = Einheitsoperator) fiir jedes pC G mit Dy =0 eine Lésung @ besitzt, und die 
Gesamtheit der Lésungen dieser Gleichung habe die Gestalt 


Gis Oy tw 


mit beliebiger Konstanten c und einem © in & abbildenden beschrankten linearen 


Operator R, fiir welchen 
IRy=0 bei ®yp=0 


iSt 7s 
Ferner gelte 
PD yo + 0 (1.3) 
os @(A,— po l)p=0 fir pew. (1.4) 
Im Spezialfall I setzen wir 
Lp = Oy = (9, Yo) (1.5) 
und 
$f 4 
Ry= f Toe ths (1.6) 
|A—po| Sd 


mit der Spektralschar FE, von Ay. Diese GréBen haben die im allgemeinen Fall 
geforderten Eigenschaften. 


1.2. Gesucht sind Elemente gm € und reelle bzw. komplexe Zahlen 4, welche 
die gestérte Gleichung 


Ayp —Ay = Bid, 9) (1.7) 
erfiillen und der ,,Normierungsbedingung“ 
Vo (1.8) 
geniigen. B(A, vy) habe die Form 
B(A, gy) = C(A, p, Aop, Ary, ---, App) = C (A, g, Aig) *. (1.9) 
Dabei bedeuten A,,A,,...,.A4 » lineare Operatoren mit dem Definitionsbereich 


und Wertebereichen in © und C eine fiir reelle bzw. komplexe 4 und g€ YW definierte 
Funktion ihrer Variablen mit Werten aus © ®. 

Die Operatoren A; ()=1,2,..., ) sollen vom Einheitsoperator J und A, in 
folgendem Sinne majorisiert werden. Es gebe fiir s,t=0 erklarte stetige reell- 


* Das Funktional / ist weitgehend frei wahlbar, R nicht eindeutig bestimmt. 
® und & kann man praktisch z. B. oft erhalten, indem man versucht, (1.2) bei zunachst 
beliebigem y zu lésen. ©@ ist Eigenelement des adjungierten Operators Aj, R bei 
Differentialgleichungen z.B. ein Integraloperator mit einer verallgemeinerten Green- 
schen Funktion als Kern (siehe das Beispiel in § 11). 

* Gleichartige Variable werden in dieser Weise mit einem Index i zusammen- 
gefaBt, 7 durchlaufe dabei jeweils alle vorher noch nicht aufgefiihrten Werte. Zum 
Beispiel ist C(A, y, Ay y, 4; y) =C(A, 9, A; 9). 

®* Der Fall, da C nur auf einer Teilmenge dieses Definitionsbereiches erklart ist, 
wird in Zusatz 2 in § 6 behandelt. 
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wertige Funktionen Q,(s, ¢) ({=1, 2,..., p) mit den Eigenschaften: 


0SQ,(s, 4) SQ,(s’, #’) tite == 6.6 FO Sr 
Q; (as, at) =a Q,(s, t) lta Orns; t 2 0, (1.10) 
| (lel A4ogll) fir per. 


Ferner moége es eine fiir alle nichtnegativen Werte der reellen Variablen 
0°, ot, oe, ...,0°*? erklarte stetige reellwertige Funktion C (0°, 01, ..., o°+®) = C (0') 
— eine Majorante fiir den Operator C — geben derart, daB gilt: 


C(0,0,...,0) =0, 

OSC +o) —C(oe) SC(H +n) -—C(e) fir 0S SE, OSo'<4f, 

ICA, p Ary) —C(wy, A:¥)| (1.11) 

<C(|A—w|+|u—m). lA: —y)||+ ]]4:¢—@)|]) — 
—C(|u— pol, , ||Aiy—o)||) fir gp per 


mit einem Bezugselement w€X, fiir das im folgenden das Ausgangselement 
des behandelten Iterationsverfahrens verwendet wird. 


§2. Die Raume ¥i und % 


Es sei ft die Gesamtheit der Kompositionen® u = ae) mit rellem bzw. 
Pu {2 
komplexem A und gwE% und YX die Menge der reellen (f+ 3)-dimensionalen 
Vektoren 0 = (0) mit den Koordinaten 9' (i =0,1, 2,...,-+2). Jedem w€ Kt wird 
ein Vektor )A| 


all Al 
»[u] =| |l4eel -( lel Jen (p20)? 
I4rell F \ Ace 
‘|p| 
als (verallgemeinerte) Norm zugeordnet. Mit 
Je] = Max|¢'| 
ee =P ell (241) 


sind in 93¢ und % numerische Normen definiert. 

MN ist damit ein Banachscher Raum. Wir beweisen, daB auch der Raum Vt 
beziiglich der in ihm definierten numerischen Norm (2.1) vollstandig ist. 

Dazu werde zunichst gezeigt, daB der Operator A, abgeschlossen ist®. Es 
sei y, (n=1, 2,3,...) eine konvergente Folge aus Y, fiir welche auch die Folge 


6 Bei dem Element u lassen wir die Indizes an A und o haufig weg. 

7 y[w] habe also als Koordinaten mindestens |A|, |||] und || A, ¢@l| und au8erdem 
alle ||A;q|| (¢=1, 2,...), welche in (1.9) wirklich vorkommen. Es k6énnen weitere 
||. A; || hinzugenommen werden (wie z.B. in Nr. 9.2). 

8 Im Spezialfall I braucht dies nicht bewiesen zu werden, da jeder selbstadjungierte 
Operator abgeschlossen ist. 


30* 
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Ao@, konvergiert. Wir setzen 
g=lmg,, @=lmg, 
mit 
8, = (Ao — Mol) On = (Ao — Mol) (Pn — 2 Gn Yo) - 

Wegen (4.1) und (1.4) ist 

L(y —2 n° Yo) =O und P#8,=0 
und damit 

Pn — LQn* Yo= Rd,,. 
In der Grenze n—> co erhalt man 
p—lo-yw=RI und OF=0, 
also ist y € Yt und 
(Ay — Mol) p = (Ao — MoZ) (VP —2P- YH) = ie 
so daB sich 
lim Ay, = lim (Ay — Mol) &, + Uolim y, = 0 + Mo” = Ao 


ergibt. Der Operator A, ist also abgeschlossen. 


Es sei nun u, (m=1, 2,3,...) eine in sich konvergente Folge aus i, d.h. es 
gelte lim |. — Un | = 0. Das bedeutet 
foe) 


[Ap Agi O10 |[Pa— Gall > 0, “Appel 04 2a pee Oe 


Daraus folgt wegen der Abgeschlossenheit von A,, daB die Folge , gegen ein 
Element » € 2% und die Folge Aygy, gegen Ayp konvergiert. Aus (1.10) ergibt sich 
dann auch die Konvergenz der Folgen A,;, gegen A;y(t=1,2,...,p). Die 
Komposition # aus A=lim/, und dem Element ist daher Grenzelement der 
Folge wu, im Sinne der Norm (2.4). Auch & ist also ein Banachscher Raum. 


yi und S haben damit die in [16], § 6 geforderten Eigenschaften. 


§ 3. Herleiten des Iterationsverfahrens 
3.1. Es seien Gy, G, und G die durch 
Gyu=A,yp—Ay, Gu=—Bid,yv), G=G,+G, (3.1) 
auf der Menge D der Elemente 
u=(*Jem mit fq = 4 
f 
definierten Operatoren. Die Aufgabe 
Gu=0 


ist dann dem Ausgangsproblem (1.7), (4.8) aquivalent. 
Der Operator Gy besitzt fiir jedes w€D die auf der abgeschlossenen Menge @ 
der Elemente 


n=(*)ex mit 1g,=0 
Pp 
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definierte Ableitun ; 
8 Gow = Ao Gn — Au Pn — An Qu- (3.2) 


Wir betrachten Iterationsverfahren der Form 
Go Uy, +Go (x4) (Mn+ 1—Un) +G, U,=G Uy t+Go(1,)(Un+1 —Uy) =0 (n=0, 1 ) 2, nh 2) : (3.3) 


Der Unterschied zum gewoéhnlichen Newtonschen Verfahren besteht darin, daB 
nur der Operator Gy (bis zum linearen Glied) nach Tay Lor entwickelt wird, 
nicht der Gesamtoperator G, und auBerdem die Ableitung an der Stelle x,, statt w,, 
genommen wird. Die x, seien dabei aus der Menge % der Elemente 


x= (M0) mit “19,45 °Ogs=e0, 
Px! 


und zwar wahlen wir y,,=¢,,°. Ausfiihrlich lautet dann (3.3) mit = a 


(An to2) Yio = (An+1 — Lo) PY, + BLA, Pn)» 
aA (EOP 2, ack (3.4) 
PoCU, lp=1, Opy+-O. 


Wendet man auf die erste dieser Gleichungen das Funktional ® an, so ergibt 
sich die Formel fiir 4,,,.,° 


ints = Ho oP B (Ay, $y) essa Ua Fe area) be (3.5) 


Nach dieser verhaltnismaBig einfachen Vorschrift (3.4), (3.5) werden die Nahe- 
rungen d,,, y, praktisch berechnet. Im Falle ® =1 ist in (3.5) ®p,=1 zu setzen 


3.2. Fiir x € & besitzt die Gleichung 
Cox) # = (Ao — Mol) Gi — An Pe = Y (3.6) 
bei beliebigem w€ © eine eindeutige Lésung hE ¥. Zunachst erhalt man 


indem man @ auf (3.6) anwendet. Wegen ®(y-+ A, y,) =0 ergibt sich dann 


Py | 
Aa Rly oo, w). 
[Goi] ist also fiir «€ § auf S durch 
os EY 
, ts D 
[Goalty = ‘ 
Dy 
definiert. Das Verfahren (3.3) laBt sich damit auch in der Form (0.4) 
Uy +4 = Uy, — [Geter Gu, (m = 0,1, 2,.. :) (3.7) 


schreiben. 


9 Unter den Voraussetzungen des Satzes in § 6 liegen die Elemente 9,,= ¢, in %, 
d.h. es gilt fiir sie auch Op, + 0. 
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§ 4. Ermittelung von Majoranten 


4.1. Bei der Untersuchung von Iterationsverfahren der Form (0.4) in [16] 
wird vorausgesetzt, daB die durch 


Hy = [Gow))*Go, Hy = [Go(w)]* G1 


mit einem geeigneten Bezugselement w definierten Operatoren Hy und H, in 
bestimmten Sinne durch im Normenraum i erklarte Operatoren H, urd H, 


majorisiert werden. Solche Majoranten werden im folgenden fiir w= %)= G a 
ermittelt. Po 

a; (i=1, 2,...,p+2) und f;(¢=0, 1, 2, 3) seien Konstanten, mit denen die 
folgenden Ungleichungen gelten: 


ean 
A Ry Sasellyll @€=0,1,2,... 


oe,| SBollyll ee | 
yweS mit ly=0O, 
| 
| 


fir ypeCS mit Gy=0, 


lv — 2° ool Al 
Py 
ay | =F 
on fiir alle ypcS. (4.1) 
|v — oe | < Pally 


Zum Beispiel kann man allgemein 
t= 1+ |Mo| o%, hire = Qj (4, %) (¢ = 1, 2, 85D) 
verwenden, denn es ist fiir pe S mit PDy=0 
AyRy = (Ay— mol) Ry tmRy=yptmoRy, 


|4oR yl S (1 + | Ho] %) Ilyll. 
und aus (1.10) folgt damit 


IA Ry SQ(IR yl, |l4oRyll) SCs.) lvl = 1,2,-..,)- 


Im Falle /= @ (also u.a. im Spezialfall I) kann man, wie unmittelbar zu sehen, 
auBerdem 
Bo == 1) ) py —— 1 


setzen, und im Spezialfall I ferner 


also 


ses y D Bs=1, Ps=1+||Po— voll 
benutzen. Denn im Spezialfall I gilt wegen ss 
(Rl=s. 
und fiir y€ S erhalt man 


2” | = (ey. vol < lvl. 
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= |lp —(Y, Yo) Poll 


ly—= Pa a;) 
S|[v — (p, Po) voll +1. vo)| Ivo — oll S (4 + |l¥o — voll) |v] 


Im Spezialfall I bereitet die Ermittelung der Konstanten «, und f; also keinerlei 
Schwierigkeiten. 
(Die Forderungen (4.1) lassen sich abschwachen. Die Ungleichungen brauchen 
nur fiir Elemente y der Form 
y=B(d,~)—B(u,y) (A, reell bzw. komplex; g,y7 EX mit lp=lyn=1) (4.2) 
zu gelten.) 
Mit diesen Konstanten «,, 6; schatzt man nun ab: 
r[(Gie)2¥] Sllvll-v fir yeS mit ly =o, 
v[(Gow) *¥] S|lpl|-%1 fir pes, 


worin bg, b, die Vektoren 


Bo Bs 
ot By a Bs 
a cer a =| oa Bs 
tp +2 Br %p+2 Bs 


bedeuten. Ferner ist nach (3.2) und (1.11) mit w= gp 
| [Gow = Gow] h|| = | A, == i, | eal “is | A,| ll Pu ie Poll, 
|G. — G,v|| = ||B(A, y)— Blu, )|| 
<C ([A—p| + | up], Lg MEN eae A (p— Il + 14: — go)ll) — 
—C(|u— po, |lp— oll, [|4i@e —go)ll)- 


Insgesamt erhalt man fiir w,v€D und hex 
» [How — How) #] S Hoow—a+yte—wp past wy) »[h], 
» [Hu — Hyv) <A, (v[u — 0] + »[v — w)) — Ay 0[v —w)) (4.3) 
mit den auf = durch 
Hye=eo-%, He=C(c)-», 
erklarten Operatoren H » und H,. 
4.2. Fiir Elemente é € %t mit °=0 ist 
[Ho plo = E109 (HB) *- 9 GG = 1,23...) 
In der Menge * der Elemente 


EEN mit &=0, #20 und £8, <1 
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existiert daher fiir o € ¥ 


=e ieee 1 
[I —Ao@)to => [Holo =9+4=ag Fe Yo- 
j=l 


Die Operatoren Hund H, besitzen die in [16] (Nr. 6.1, 3.1) geforderten Eigen- 
schaften, und alle dort vorkommenden Begriffe sind erklart. 


§ 5. Vergleichsaufgabe und Vergleichsverfahren 


5.1. Zur Konvergenzuntersuchung und Fehlerabschatzung haben wir nach [16] 
(Satz und Zusatz 4) eine Vergleichsaufgabe mit dem Vergleichsoperator 


G a Go — Gy 
zu lésen, wobei Ge und G; die durch 


Goo = [I —Hogy] a r Qo Ho ot+Hoo, 


Gy G= —H 100 5 H 1@Q 
auf 9 definierten Operatoren bedeuten und hier 
| Ao — Ho 
0) 
§ =v[% — w] =0, Oo =¥[% — vw] = 0 
0) 
sowie 
| Ay — Ao| 
» [U, — Up] = ler — oll So, = (7) EN 
|4 5 (G1 — %o)| 
ist. Gesucht wird ein Element t= 0, aus Jt mit 
Grso, (5.1) 
also etwa eine Lésung der Gleichung Go =0. 
(5.1) lautet ausfiithrlich 
A+ e+ Bott BrlC(r}) —C (0)] 2°, 
py a Be 
o1 t o,f, Tt + opBalC (t!) —C(oo)] St (k= 1,2,...,p +2) ine, 


Sind 7° und t! Zahlen, mit denen 


Ate+ Art BZ, 


C(r9, ah a+ S (F#—oi))—C(@8, 0,0,..., 0)] S 2°, 


1 + o% By 9 t+ 0% Bs IE (=*, 7, Ce a (—o}))—C (6, OGe 0)) <7 


1 


gilt, wobei z die Werte 2, 3,..., #-+-2 durchlauft, so erfiillen die Zahlen 


TOT et Oy aes 04) (0 ="253 pe 2) 


die Ungleichungen (5.2). Zur Ermittelung von 7 geniigt es also, die beiden nicht- 
linearen Ungleichungen (5.3) zu lésen. 
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5.2. Mit den Elementen 


lautet das Vergleichsverfahren (0.6) 


Ont+1 = On — [Go,] 2G On (n =0,1, 2,.. .) 


koordinatenweise geschrieben 


Char = — par fol + ob + Bell (oh) —E (0d)]}, 


Bo On 
On¢1= 91 + 0181024104 + a Bs (C (04) —C (o4)], (5.4) 
Oi = a (OG 0) (B= 2,3,:-.,Pp+2) (n=0,1,%..). 


co 
Soweit die &, aus § sind, dh. BoOn<1 gilt, geniigen die Elemente w,, x, und 
0,,&, den in [16] gestellten Forderungen 


W [Uy ae: Ky | = Cn _ Ey (5 5) 
5. 


Pg ay ey ee ene (n= 0, 1, 2) 6-35 X_1 = Xo, E,=0). 


Diese Ungleichungen folgen mit Hilfe vollstandiger Induktion daraus, daB nach 
[16] (Nr. 2.4) 
» [Uy — Up—1] S Ox — Ora (5.6) 


gilt, soweit die vorkommenden Elemente mit (5.5) geniigenden Elementen x,, 
und &,€% erklirt sind. Man vgl. auch das Beispiel in [6] (§ 4), wo ahnliche 
Elemente x, und &, verwendet werden. 

Besitzen die Ungleichungen (5.3) eine Lésung 1°, t! mit Byt!<1, 1°9= 0}, 
ti=o}, so folgt aus dem Zusatz 4 in [16], daB die Iteration (5.4) unbeschrankt 
durchfithrbar ist und e,<t (unter anderem also By0;<fot'<1) gilt. Da die 
Vektorfolge o,, auBerdem monoton zunimmt (vgl. (5.6)), konvergiert sie dann 
gegen einen Vektor o*<t. Dieser list die Gleichung Gg =0, was man einsieht, 
indem man in (5.4) den Grenziibergang  — oo durchfiihrt. Die dann entstehenden 
Gleichungen sind mit Go*=0 identisch (vgl. (5.2)). 

Existiert also eine solche Lésung der Ungleichungen (5.3), so sind alle Voraus- 
setzungen des Satzes 4 in [16] erfiillt, und es ergibt sich das im nachsten Ab- 
schnitt formulierte Ergebnis. 


5.3. Am genauesten werden die Fehlerschranken, wenn man o,=v|[u,— Uo] 
setzt oder v[u,;— | im Einzelfall moglichst gut abschatzt. Jedoch wollen wir 
diese Norm auch allgemein abschatzen, miissen dabei aber einen u.U. erheb- 
lichen Genauigkeitsverlust in Kauf nehmen. Eine solche Abschatzung fiir 
y[u,—U,| wird z.B. bendtigt, wenn man q, gar nicht berechnet hat. 
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Man erhalt aus (3.4): 


ates | bo — Ao + Ht mit ,=— éx D B(Ag, Yo); 
Yo— Pot Rta Go + B (Ao, Po) ] ibe 
oe [eo — ol + |e 
v[%,— |S ll Po — oll + 9 mit #=%4|1Po+ B(Ao, Poll. (5-7) 
A; (Po — Yo)| AWE a UY 


Auf der rechten Seite dieser Ungleichung kommen nur GréBen vor, welche 
man ohne Kenntnis von q,, also ohne eine Gleichung lésen zu miissen, berechnen 
kann. 

§ 6. Konvergenzaussagen und Fehlerabschatzungen (Ergebnis) 


6.1. Wir fassen die aus den bisherigen Untersuchungen und [/6] (Satz und 
Zusatz 4) folgenden Ergebnisse iiber das Iterationsverfahren (3.4), (3.5) 


1 
Ana = bo — toe Pn)» 


(Ag — Mo!) Pnt1i = (An+1 — Mo) Pn + B(An» Pn) > lOn+1 = 1, (6.1) 
Moet mit lgy=1, Dy+=O PO, Toe 
(Lp =®ey =(¥, Yo) im Spezialfall J) zur Lésung der Aufgabe (1.7), (1.8) 
Ayp—Ap=Bii,), lo=t (6.2) 
zusammen. Die in §1 genannten Voraussetzungen seien dabei erfiillt. 


Zur Fehlerabschatzung hat man folgende Rechnungen durchzufiihren: 
1. Man schatze 
a) die Operatoren A,(t=1, 2,..., ) in der Form (1.10) und 
b) den Operator C in der Form (4.11) mit @ = @p ab. 
2. Man ermittele mit Hilfe der Tabellen 1 und 2 die Konstanten «;(7=1, 2, ..., 
p+2), B;(¢=0, 1, 2,3) sowie die Schranken 0§ und of (¢=0,1, 2,..., p+ 2). 
3. Man stelle fest, ob die Ungleichungen (6.3) eine Lésung besitzen. 


Satz. Gibt es Zahlen 7°, t! mit P=S>o}, 1201), Byt!<1, welche die Unglei- 
chungen 


O+et Bort Bo ie (e, 0% = (7 —o}))— C (08, 0, 0, a) SDE) 
1 


, | : (6.3) 
O1 + o B, 7°t+ a, Bs ie ee th of 4 = (z} o1)| CXOny OnCan ss 0)| = 

1 
erfiillen (wobet 1 die Werte 2,3,...,p+2 durchliuft), so konvergieren™ die dann 


durch die Iterationsvorschrift definerten Folgen 2,, p, gegen eine Losung A*, y* der 


© Es geniigt, statt to} die Ungleichungen t!=>o!— “ of (¢=1,2,...,p+2) 
nachzupriifen. t!= oj folgt dann aus (6.3). re 

“Wenn wir im folgenden von den ,,Gleichungen (6.3)‘‘ sprechen, so sind damit 
die Beziehungen (6.3) mit = statt < gemeint. 

“ Die Konvergenz von Elementenfolgen ist dabei im Sinne der jeweils verwen- 
deten Norm zu verstehen, lim y,= y* bedeutet also lim ||p*— gy, || =. 
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gegebenen Aufgabe, und es gelten die Fehlerabschatzungen 


ee ee eet (6.4) 
Una allgemelner 


[AF —A,|Se—eh, llo*— gill St—e,  (w=0,1,2,...) (65) 
mit den durch (5.4) definierten Zahlen 0°, 0}. 


Zusatz 1. Unter den Voraussetzungen des Satzes besitzen die Folgen A; @, die 
Grenzelemente A, p*, und es gilt 


|4s(o* — o)|| SA  — @) (mn =0,1,2,...;7=0,1,2,...,~). (6.6) 


Zusatz 2. Es sez D, ein offenes Intervall der h-Achse bzw. ein offenes konvexes 
Gebtet der A-Ebene, in welchem iy liegt, und D, eine offene konvexe Teilmenge von 
Elementen p CX mit lp=1, welche wo und qo enthdlt. Ist B(A, py) =C (A, 9, 4; ¢) 
nur fiir Zahlen KCD, und Elemente p€D, erklart und die in (1.11) enthaltene Un- 
gleichung fiir C in dieser Teilmenge erfiillt, so gelten die Aussagen des obigen Satzes, 
falls D, alle Zahlen 2 mit |A—A,|S1°—oj— 09 und D, alle Elemente pEX mit 
lp=1 und ||p—|| Sc—ot, ||4:(@—)|| S$? (Co) @=0,1,2,..., 9) 
enthdilt. : 


Die Richtigkeit dieses Zusatzes 2 ergibt sich, wenn man in [16] als D die 
Menge der Elemente w€ ® mit 1,€9,, 9,€D, und lp, =1 verwendet. 


Bemerkungen. Ist A, beschrankt, so kann man Q;(s, ¢) = || A;|| s verwenden. 
Fiir nichtbeschrankte Storoperatoren bei gewohnlichen Differentialgleichungen 
sind Abschatzungsfunktionen Q; in [14] angegeben (s. auch § 9). — Majoranten G 
werden im folgenden fiir verschiedene Falle hergeleitet. — Im Spezialfall I er- 
fordert die Ermittelung der Konstanten «, und 8, nur sehr geringe Rechenarbeit. — 
Fiir den Spezialfall I und ein Stérungsglied der Form B(A, py) =—A,y werden 
die Gleichungen (6.3) in § 7 gelést. Ist A, dabei hermetisch, so lassen sich die 
Fehlerschranken noch verbessern, wie in 7.3 gezeigt wird. 


6.2. Es sei noch bemerkt, daB man die in dem Satz enthaltene Existenzaus- 
sage sowie die Fehlerabschatzung (6.4) fiir die erste Naherung auch auf einfache- 
rem Wege bekommen kann. Das untersuchte Iterationsverfahren stimmt fiir 
n=0 mit dem Verfahren 


An+1 = Lo ae Di, Lo) (Pn — Po) a aa OB (An , Pn) , 


(Ap— Mol) Gnti= (An+1— Mo) Po + (A,—!Mo) (Pn— Po) + Blan» Pn)> $Pn41=1 (6.7) 


MEU mit lp=1, Ogy+0 (7 =O eer) 


iiberein, d.h. die jeweiligen ersten Naherungen sind gleich. Dieses letztere Ver- 
fahren kann man aber in der einfacheren Form 


tn Lu, (te OAD.) 
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mit dem durch 


1 
DBA, 
rade) 
Tu=| pot R\(A— Mo) (P — Po) ee D(A — fo) (P — Po) Po + 
+ BUL 9) — gq PBL) 


Mo — aa D(A — bo) (YP — Po) — © 


auf D definierten Operator schreiben. Darauf laBt sich dann unmittelbar der 
Satz 1 von [16] anwenden, welcher solche Verfahren als Spezialfall mit behandelt. 


Durch im wesentlichen die gleichen Abschatzungen wie im § 4 beweist man 
eine Ungleichung 
y[Tu, Tv] SF lo[u—v],»[u—w],r[v—w]| mit w= ka ee 
0 


der in [16], § 1 geforderten Form, und Satz und Zusatz 1 aus [16] ergeben dann, 
daB fiir das Verfahren (6.7) der oben formulierte Satz gilt, wenn die @,, jetzt die 
Naherungen des zu (6.7) gehérenden Vergleichsverfahrens 


On+1 = 97 55 00 ae 0}, On Do se [C (oi) = C (oi) ] Dy (n = 0,1, 2, Bec) 
bedeuten. 


Andererseits bekommt man diese Ergebnisse fiir das (im allgemeinen lang- 
samer konvergierende) Verfahren (6.7) auch durch geringe Variation unserer bis- 
herigen Uberlegungen. Man braucht nur x,=x,) und &,=0 (n=O, 1, 2,...) zu 
setzen. 


§ 7. Spezialfalle bei linearen Eigenwertaufgaben 


Wir untersuchen einige Spezialfalle der allgemeinen Aufgabenstellung genauer. 
es eal sito 
BA, 9p) =— Ap + Aap 


geht die Gleichung (1.7) tiber in 
(Ay + A))~p=AUL + A,) 9, 
und man schatzt dafiir ab: 
|| B(A, @) — Blu, y)l| S|]41 @ — Aryl] + [ol | 4ee — Aeyll +14 — | || Ao gol] + 
a A el 43) 5) eg |) 
x [(l|42¢ — Aayll + || dep — Ae qoll) — |] 4ey — Aa@oll] + 
t= [([4—p] + |¢ — mol) — | — mol | x 
x[(l|42@ — Asyll + ||42v — Aegoll) + ll 4ey — 4s goll]- 


Also gilt (4.44) mit 
C(o") = e + | ol o* + 0° || Ao ol] + 0° of. (7.1) 


Fiir 4,=0 erhalt man die besonders einfache Majorante 


C/o. (7.2) 


Eigenwert- und Verzweigungsaufgaben 453 


7.2. Es seien nun die Voraussetzungen des Spezialfalles I erfiillt, und auBer- 
dem habe B die Form 
Bid, 9) =— Ax. 


Damit haben wir den in [13] behandelten Spezialfall vor uns. Hierfiir lauten 
die Gleichungen (6.3) 


+a [c—b+ y] =~, 
(7.3) 
DEY ep Oe) ty, 
wenn ; : 
und 4 
arahers tos Ae 
gesetzt wird!4, Man kann sie leicht lésen und erhalt das 
Ergebnis **. Gilt (im Spezialfall I bei B(A, y) = — 4,9) 
1—(aB+actab+a\)=0 
und : (7.5) 


r=0 mit r=|1— («B+ a(c —b) + a)]?—4a(b+aB(c — By), 


so konvergieren die dann durch (6.1) definierten Folgen i, , p, gegen eine Lésung 
A*, p* der gegebenen Aufgabe (6.2), und es gilt die Fehlerabschatzung 


[ae A lsa@ a), (7.6) 
let*-gillSy—b, ||Ao(e*—e)ll S 44+] 01) (v—2), |] 41 (9*-q)|| S da (y—2) 


mit a 


t=abtelcab-bay, y= (Au Dee edly 


Setzen wir nun speziell Ap= Wy, Po= Yo. SO ist A, mit der durch Potenzreihen- 
entwicklung berechneten ersten Naherung “y+, fiir den Eigenwert identisch, 
so daB wir die entsprechenden Fehlerabschatzungen vergleichen kénnen. Es wird 
nach (5.7) 


[Ai — Ao| = | “44 = || Ar voll = 0: (||o 


’ 


Aol) 7 Q1(1, | ol) =ydsad, 
(7.7) 


lla: — Moll SO | 1 : 1 
G=—_J\|A 2_ 14, |2?< —||A Sy IO, 
mit 1 
ba = (Ar Yo. Yo)» y=alz. elke a= (7, 1 eet) 


Verwendet man diese Zahl«, B=1 und als of, of und oj in (7.4) die Schranken 


1 
Ha], a || Ai pol] und « || 41 yo 
Schranken sowie 09 = 0, so erhalt man eine Fehlerabschatzung (7.6) fiir A,, welche 
in [13] auf anderem Wege bewiesen wurde. Die entsprechende Fehlerabschatzung 


mit der genaueren GréBe # leitete F. W. SCHAFKE in [12] auf wieder andere 
Weise her. 


| oder irgendwelche der in (7.7) stehenden gréberen 


14 Den trivialen Fall « =0 lassen wir dabei auBer acht. 
15 Diese Verbesserung lat sich auch durch geringe Variation des Beweises in 
[13] erreichen. 
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7.3. Fiir den in 7.2 behandelten Fall lassen sich die Abschatzungen verbessern, 
wenn A, Hermetisch ist. Es ist dann 
(Ay (y —), Yo) 


— H,v=— [Go| 2 (BA, g) — Blu, -( 
Face Haver Gol Ae Phe att als AB ne i ea 
und 
(Alp —), vol =p — Ai — 12) Yo)| S [lo — vl I(41 — 4D) voll = 44 lle — vl 
(9 wie oben), so daB die Ungleichung (4.3) jetzt mit 


md 0 o AS 0° 
fo ( Js [stat fhe=(. 50) 


gilt (C = 9% nach (7.2)). Statt der Gleichungen (7.3) erhalt man 
at+dy=x, b+xy4+ap[e—b+y]=y 


und damit das 


Ergebnis (fiir den Fall 9==0). Ist (im Spezialfall I bei B(A, py) =— Ay) Ay 
Hermitesch, so gilt das in 7.2 formulierte Ergebnis auch, wenn man die Konver- 
genzbedingung (7.5) durch 


1—(«B+20b+a)=0 und rZ=0 mit r=[A— («B+ a)]}?—40[6+ a8 (c—d)] 


ersetzt und die GroBen Se 
K=atOyY, Y= iit ores ts 


20 


verwendet 1°, 


Auch in diesem Ergebnis sind wieder die fiir Hermiteschen Operator A, in 
[13] bewiesene sowie eine von F. W. ScHAFKE [12] hergeleitete Fehlerabschat- 
zung fiir A, enthalten. 


§ 8. Eigenwertaufgaben bei Matrizen 


8.1. Zu lésen sei die Eigenwertaufgabe 
(My +My) gp = A(No + M4) g, (8.1) 


worin g einen m-dimensionalen Vektor bezeichnet und M,, M,, Nj und N, mx m- 
Matrizen bedeuten. Die Matrizen My und N, seien Hermitesch (My= Mo, Ny= No 
und N, sei auBerdem positiv definit (p’ Nyy > 0 fiir m+ 0). 
Die ,,ungestérte Aufgabe“ 
My =ANY (8.2) 


besitze einen einfachen Eigenwert jy. wo sei ein zugehériger durch 
Po No Po = 1 
normierter Eigenvektor. Gesucht sind Lésungen A, m der Aufgabe (8.1) mit 
Po Nop = 1. (8.3) 


*° Die sich hiermit ergebenden Fehlerabschatzungen sind im allgemeinen genauer 
(oder jedenfalls nicht ungenauer) als die aus dem Ergebnis in Nr. 7.2 folgenden. 
Sicher ist dies z.B. der Fall fiir c—b 20. Oft wird c=b=0 sein, vgl. (7.7) 
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Jeder weitere Eigenwert des Problems (8.2) habe von jg mindestens den 
Abstand d. Mit x bezeichnen wir den kleinsten Eigenwert der Matrix N,, mit 
|M,| und |N,| nichtnegative Konstanten, fiir welche 


(M,¢) Mos IM,"?o >, (Moy Mes IN|? o' y (8.4) 
gilt. Sind g;, die Elemente der Matrix M,, so kann man z.B. 
|M,| = |Z lan? (8.5) 
i,k 
oder 
Wap Seo 2 Max > lanl, S=Max >| gus! 


verwenden (entsprechend bei N),). 

Wendet man auf diese Aufgabe den Satz aus § 6 an, so erhdlt man das am 
SchluB der Nr. 8.2 formulierte unmittelbar anwendbare Ergebnis. Ein darin 
enthaltenes spezielleres Ergebnis ist auBerdem in Nr. 8.3 genannt. 


8.2. Dieses Problem laBt sich als Aufgabe der Form (1.7), (1.8) deuten, bei 
welcher die Voraussetzungen des Spezialfalles I erfiillt sind, und wir kénnen dann 
den Satz aus §6 anwenden. © sei der Hilbertsche Raum der Vektoren p mit 
dem inneren Produkt 


(vv) =~ Ny 


IIe =VO Nog. (8.6) 


und der Norm 


Die Aufgabe 
Ayp—Avg=BiA,¢), (Yo) =1 


mit 
B(4,¢) = —Ay~p +4429, Ay, = Noe My; Ayo M;, Ase=NesN, 


ist dem zu lésenden Problem (8.1), (8.3) aquivalent. A, ist dabei selbstadjungiert. 
Das zugehorige Iterationsverfahren (6.1) 1a4Bt sich in der Form 


Antt = Ho + Po M, Pn — ‘1, DN. 1Pn»> 
(My — Mo No) nti = (Anta — Ho) No Pn — My On + An i Pn (8.7) 
wo NoG=— 1, BoNo Prai—1 Op oe) 


schreiben, wobei /p = ®y =(¢q, Wo) gesetzt ist. 
Wir schatzen B(A, y) in der erforderlichen Weise ab. Man erhalt 


lA. 9P=Oh 9) No Me S— (9 MoS —|MP9' ee ltEle Re 


(8.8) 
|4e@|P2= (Me) No” Nes 9) Np s— [M2 9’ gS IM l2o Noe, 
also i 
M 
WA, al oll, 4eglls Hol (8.9) 
und damit die Zahlen 
M. N, 
wage antl, alt well. sn 
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Als € laBt sich nach 7.1 die in (7.1) angegebene Funktion 
C (0°, ot... of) = 08 + [Ho] o* + 0° || 42 Goll + 0° 0" (8.41) 


oder wegen (8.8) auch 


~ 


N, - 
E (0%, hs. 08) = 0° + [uo] ot + 0° ML I poll + ore! (||Pol|=1 fiir po=Yo) (8.12) 


verwenden. 


Ergebnis. Fir das Iterationsverfahren (8.7) zur Losung der Aufgabe (8.1), (8.3) 
gilt (unter den in 8.1 genannten Voraussetzungen) der Satz aus $6 mit p=2, wenn 
man darin 


als ||p|| die Norm (8.6), 

als C die Funktion (8.12) (oder (8.11)), 

als a; die Zahlen (8.10), 

als B; die in Tabelle I fiir den Spezralfall I genannten Konstanten, 

als \oos a, (t=0,1,3,4) die in Tabelle 2, Spalte 2 erklarten und unter Be- 
nutzung von (8.8) zu berechnenden Schranken verwendet. 


Die Gleichungen (6.3) sind hier leicht zu lésen, sie fiihren auf quadratische 
Gleichungen fiir 7° und 71. 


8.3. Wir leiten aus diesem Ergebnis noch eine einfache Fehlerschranke fiir 
|A*—uo| her, welche nur die oben erklaérten Konstanten d, x, |M,|,|N,| enthalt. 
Sie kann etwa dazu dienen, den Einflu8 von Abrundungsfehlern bei den Matrizen 
M, und N, auf den Eigenwert wy der Aufgabe (8.2) abzuschatzen. 


Es ist 


[(B (Mo, Po)» Po)| = || B (uo, Yo) ) || B (uo, Yo) — (B (Uo Yo), Yo) Poll = I| B (Ho Po)| 


und 


Bo, voll S |Asvoll + [Hol | Aavoll S 9-4 


g=— > (IMi| + [ol |N4l). 


mit 
Aus (5.7) folgt damit fiir Ag= py, Yo=Yo 


Yt — y's oy h Init eles Geld | Men 


(q+sx)(Aty)=%*, (q+sx)(Aty)=y—xy bei pe ae y=, SEAL 


Gilt s+2)/¢S1, so existiert eine Lésung x, y mit 


_i-—s Gs)? 2q 
ve 2 / 4 SES (8.13) 


Der Satz aus § 6 liefert dann das 
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Ergebnis. Gilt 


d|N,| +2 xd (|M,| +] mo] |M|) Sx, (8.14) 
so besitzt die Aufgabe (8.1) einen Eigenwert A* mit 
M N, 
|A* — 11| <2 tel i! (genauer: | A* — uo| dx). (8.15) 


Ist Ny die Einheitsmatrix und N,=0, so gehen diese Formeln (8.14), (8.45) 
uber in 
4|M,|Sd@ bzw. |A*—po| S2|M]| 


d = Vaal M 
[genauer: | A*—o| S 5 (! \1 4! ze) Mal (14 Sl), 


\ 


Fiir die Eigenwertaufgabe 
(M+ M,) p=Ap 


mit normalen Matrizen M, und M, bewiesen A. J. HOFFMANN u. H.W. WIELANDT 
[4], daB allgemein (also ohne Zusatzbedingung) 


Da — 6? SIP (8.16) 
gilt, wenn yi die Eigenwerte der Matrix M, und /*? diejenigen der Matrix M,+M, 
bedeuten und man diese Eigenwerte in geeigneter Weise numeriert. ||| 14,||| ist 
dabei die in (8.5) genannte Konstante |M,|. Wie Horrmann und WIELANDT 
bemerken, braucht eine solche Ungleichung bei nichtnormaler Matrix M, nicht 
zu gelten. Aus (8.16) folgert man fiir den einzelnen Eigenwert 


| A* — po] SII] M4Ill- 


Zum Vergleich sei bemerkt, daB bei normaler Matrix M, mit den Eigenwerten 


x, dieser Matrix ego 
ieetl=|/ Shor 


gilt und die kleinstmégliche Konstante |M,| in (8.4) 
|M,| = Max | x,| 
ist. 
§ 9. Gewohnliche Differentialgleichungen 
9.1. Als Beispiel behandeln wir die Eigenwertaufgabe 


— p+ M[y] =Age+ P(A) N[g], (9.1) 
ap(0)—Be'(0)=0, yet) + det) =0. (9.2) 
Dabei bedeuten M und WN lineare Differentialausdriicke 
M [¢] = b2(*) y+ Pi(*) y+ Po(*) @, 
N [gy] = N(x) g + Go(*) P 
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mit stetigen Funktionen #;(x) und q;(x) und 
ps(x)<1 (0S%3S1), (9.3) 
P(4) ein Polynom in A und «, B, y, 6 Konstanten’” mit 
oe, B 7, OZOl™ cP BS 05. ty Oo 0: 


c,; und d; seien Zahlen, fiir welche (im Intervall [0, 1]) 


|p:(x)| Se; (1 = 0,1, 2), |ai(*)| $4; (1 = 0, 1) 
gilt. 
Die ,,ungestorte“ Aufgabe 
—g"=Agp, a«(0)—Be'(0)=0, yep(t)+dy'(1)=0 (9.4) 


besitze A=jp als einfachen Eigenwert, y=w (x) sei eine zugehdrige, durch 


1 
dl poey ee nd 


normierte Eigenfunktion. Jeder weitere Eigenwert der Aufgabe (9.4) habe von 
4) mindestens den Abstand d. 


Gesucht wird eine Lésung A, w(x) der Aufgabe (9.1), (9.2) mit stetigem q’’(x) 
und 


1 
[ 9) vols) dx = 1. (9.5) 
Das Ergebnis fiir dieses Problem ist am SchluB8 der Nr. 9.2 in unmittelbar anwend- 


barer Gestalt formuliert. 


9.2. Wir schreiben die Aufgabe als St6rungsproblem im Hilbertschen Raum 
GS =X,/0, 1] mit dem inneren Produkt 


(, ) =F play ln) dx 


und der Norm 


lol =|/flooneas. 0.) 


% bedeute die Menge der Funktionen q(x) mit stetiger zweiter Ableitung, 
welche die Randbedingungen (9.2) erfiillen. Ay, A, und A, seien die auf Y& durch 


Ayy=—g", A\p=M[q], 4.9=N[q] 


definierten Operatoren. Mit 


B(A, ») = — A, + P(A) Arp =C(A, Ay gy, 429) (9.7) 
1aBt sich das gegebene Problem in der Form 
App —Ap=BiA,g), (~. Yo) =1 (9.8) 


™ In § 7 wurden die Buchstaben «, B, y bereits in anderer Bedeutung verwendet. 
Eine Verwechselung ist kaum zu befiirchten. 
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schreiben. Die Tatsache, daB A, nicht selbstadjungiert ist im Sinne der in der 

Funktionalanalysis iiblichen Definition!’, erfordert im folgenden zusatzliche 

Uberlegungen. Ohne Riicksicht darauf soll zunachst B abgeschatzt werden. 
Nach [14] ist 


Ar ells (ell (409! 


iE lA2¢llS Q2(llell, | 


0,(s,t) = est c, Vstt cos, - On(s, 2 =a, Vst+ as. (9.10) 


Ay ¢||) (9.9) 


mit 


Ferner gilt fiir die durch (9.7) definierte Funktion C eine Ungleichung der in 
(1.11) geforderten Art mit wm = q, und 


C(Q°, ot, ..., 04) = @ + Plo? + | Mol) (o* + ||N Goll) — P(|Hol) IN voll (9-41) 


Diese Funktion C leitet man auf ahnliche Weise her wie die im Beispiel von § 11 
benutzte. 


Das (6.1) entsprechende Iterationsverfahren fiir die hier behandelte Aufgabe 
lautet mit /~ = Dy =(¢g, y) 


1 
Agee = fo Pe (M [Pn] a P(A,) N [@n]) Py AX, 


Tt Pn+1 — Mo Pati = (A, == Ho) Pn — M [Pn] ef P(A,,) N [Pn] , 


, 1 , (9.12) 
% Pnt+1 (0) Pr B @n+1 (0) rik 0, Y Pn+1 (1) ale 6 Py +1(1) a 0, 


1 1 
I Posed 1 f Onur Pods = 4 (Vir On Ah Dig. 2) 
0 0 


(¢o(x) erfiille dabei die Randbedingungen (9.2), und yo(x) sei stetig). 


Der Operator A, besitzt ein in © vollstandiges Orthonormalsystem von Eigen- 
funktionen. Er laBt sich zu einem selbstadjungierten Operator A, mit dem 
Definitionsbereich % durch AbschlieBen fortsetzen. % besteht aus allen Grenz- 
elementen gy=lim g,, von Folgen q,,€ %, fiir welche auch die Folgen Aggy, kon- 
vergieren, und A,p wird dann durch A,y=lim A @, definiert. Wegen (9.9) 
kann man auch die Operatoren A, und A, zu auf 9% erklarten Operatoren 4A,, A, 
fortsetzen, indem man 4,y=limA;@, (i=1, 2) definiert. Da der Operator A, 
die gleichen Eigenwerte wie A, besitzt und die Funktionen Q; und C stetig sind, 
sind fiir das (9.8) entsprechende Problem mit den fortgesetzten Operatoren alle 
im Spezialfall I geforderten Voraussetzungen erfiillt. Es gelten dafiir also die 
Aussagen des Satzes aus § 6, dessen Voraussetzungen im folgenden erfiillt seien. 

Die Elemente y,€% des entsprechenden Iterationsverfahrens sind sicher 
sogar aus Y%, wenn man g@€ % wahlt, vgl. (9.12). Wir zeigen, daB auch * in Nt 
liegt. 

Es sei w eine reelle Zahl, welche kein Eigenwert der ungestérten Aufgabe ist. 
G(x, €) bedeute die zum Differentialausdruck 


— Pp Ly (9.13) 


18 4, ist nicht abgeschlossen. 
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und den Randbedingungen (9.2) gehérende Greensche Funktion. Dann ist 
(A,—wI)1@ auf der Menge der stetigen Funktionen durch 


(A, —phag= fG(m8 y (&) dé 


definiert und (A,—J) 19 in der gleichen Form fiir y¢ S erklart. Daraus folgt, 
daB jede Funktion aus {{ als dem Wertebereich des Operators (A,—) + stetig 
ist und eine stetige Ableitung besitzt. g* geniigt der Gleichung (9.1), wenn man 
dort —q@” durch A,p ersetzt. Wegen (9.3) ergibt sich daher, daB (A,— 1) g* 
stetig, also aus dem Wertebereich des Operators (Ag—J)7 ist, so daB p* aus A 
sein muB. 

Man kann auch den Betrag | y*(x) —q,,(x)| abschatzen. Wir setzen A,= 
A,—I (u wie oben). Dieser Operator kommt in (9.7) (ebenso wie auch Ao) 
explizit nicht vor, jedoch kénnen wir den Zusatz1 aus §6 ohne weiteres fiir 
p =3 anwenden, also auch den Fehler || A,;(g*—¢,,)|| abschatzen. Man hat 


5? 


[45 R|] = Max OAT 1 4 dered _,, 


wobei R den Operator (1.6) bedeutet und 4 die Eigenwerte von Ay, abgesehen 
von A=), durchlauft. Damit erhalt man unter den Voraussetzungen des Satzes 


5 


wegen = =d-+ |y—| die Fehlerabschatzung 


|| 43(@* — ei S (@+ |4o— Hv) (A— el) = (w = 0,1, 2,...). 


Mit Hilfe der Schwarzschen Ungleichung folgt 


| e*(x) — py ( |=|f ews 3(P* — Pn) yae| s/f 1e0e. oyras 4a (p* — 9) 
(G(x, €) wie oben), also 


1 
19%) 960) =| FG. 8)8d8 2+ [ew —al) (8 0) 
(1: O12 pen} 0; 101). 


(9.14) 


Verwendet man in (9.14) t!= 9*1, so ergibt sich wegen o,—>0*, daB die Folge 
P, (x) gleichmaBig gegen p*(x) konvergiert. 


Ergebnis. Fiir das Iterationsverfahren (9.12) zur Lésung der Aufgabe (9.1), 


(9.2), (9.5) gilt (unter den in 9.1 genannten Voraussetzungen) der Satz aus § 6 
mit p=2, wenn man darin 


als ||q|| die Norm (9.6), 

als Qh, Qs die durch (9.10) erklarten Funktionen, 

als C die in (9.11) genannte Funktion, 

als o;,B; die in Tabelle 1 fiir den Spezialfail I aufgefiihrten Konstanten, 

als 05, 01 ((=0,1,3,4) die in Tabelle 2, Spalte 2 bei A,p=—Mgq, A.9=No 
erklarten Schranken verwendet. 
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Unter den Voraussetzungen dieses Satzes konvergiert die Folge ,(x) sogar 
gleichmipig (im Interval [0,1]) gegen p*(x), und es gilt die Fehlerabschitzung 
(9.14), worin u eine von den Eigenwerten der ungestorten Aufgabe (9.4) verschiedene 
veelle Zahl bedeutet und G(x, &) die zwm Differentialausdruck (9.13) und den Rand- 
bedingungen (9.2) gehdrige Greensche Funktion ist. 


§ 10. Zwei einfache Beispiele 


Die zwei sehr einfachen Beispiele dieses Abschnitts dienen vor allem dazu, 
einen (mehr theoretischen) Einblick in die Genauigkeit der Abschatzungen zu 
erhalten. Wir verwenden dabei die Bezeichnungen des § 8. 


10.1. Beispiel 1. Wir rechnen zum Vergleich die schon in [13] behandelte 
Matrizenaufgabe 


(M)+M,)p=A4p, Prp=1 (10.1) 
mit 
‘On RO O-—=4 
Mo=(‘5 a i=(( 0 ) 


und zwar soll die 
1 
Mo= 10, YPo= (3) 


benachbarte Lésung ermittelt werden. 


Die Tabelle 3 enthalt die mit der tiblichen Stérungsrechnung (Potenzreihen- 


entwicklung) erhaltenen Naherungen >) wu, fiir den Eigenwert (Spalte 2), die 
v=0 
exacten Fehlerbetrage (Spalte 3) sowie die mit den Ergebnissen von [13] berech- 


neten Fehlerschranken (Spalte 4)?9. 


Tabelle 3. Evgebnisse der tiblichen Storungsrechnung 
(Potenzreihenentwickelung nach ¢, vgl. [13]) 


Naherungen s Be 
n |A*—2,|= |A*—Ay|S 
Li Tis » Ly (exact) (nach Stérungstheorie [13]) 
v=0 
1 10 0,0501256... 0,0527782 
2 9,95 0,0001256... 0,002 7782 
3 9,95 0,0001256... 0,000 2782 
4 9,949875 0,000 0006... 0,000 028 2 
5 9,949875 0,000 0006... 0,000 003 2 


In Spalte 2 der Tabelle 4 stehen die Naherungen /,,, welche man mit fo, Yo 
beginnend nach der Iterationsvorschrift (8.7) erhalt, Spalte 3 enthalt die exacten 


: é M 
Fehlerbetrage. Die Fehlerschranken in Spalte 4 wurden nach (7.6) mit «= [Ml = x 


und B=1+||q—Yo|| ermittelt. Dabei wurden die jeweils abzuschatzenden 
GréBen als erste Naherung 2,, y, aufgefaBt, als gy also die vorletzte Naherung 


19 Wegen ,—0 ergibt sich hier bei Verwendung der Gr6éBe @ in (7.7) keine Ver- 
besserung. 
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verwendet und Aj=y=10 gesetzt (Ay kommt hier in der Iterationsvorschrift 
nicht vor, ist bei der Abschatzung also frei wahlbar). 


Tabelle 4. Ergebnisse bei Anwendung des Iterationsverfahens 


; |a*— ni = |a*—ay|S 
Naherungen : 

n a (exact, zugleich Schranken (nach (7.6) mit 

Hn a nach (7.6) mit 6=|g8)) B=1+ |lro—vall) 

| | 

1 10 0,050125629... 0,0527782 
2 9,95 0,000125629... 0,000 1323 
3 9,949875 0,000 000629... 0,000 000 663 
4 9,949 874374... 0,000 000003... 0,000 000 003 4 
5 9,949874 371... 


Das Iterationsverfahren liefert bessere Werte, und die zugehérigen Fehler- 
abschatzungen sind genauer. Daf die ersten der genannten Fehlerschranken 
iibereinstimmen, hangt mit der Einfachheit der Aufgabenstellung zusammen. 


Andert man in der Aufgabe M, ab zu 
A eC) ae 
so wird 


Perec 1 
— a 2 — fh Nie Be b 
A*=]/100—2?, y= ee — 100 — =) 


Bar A, os Co— oe — as und f=1 ist die Konvergenzbedingung (7.5) 


erfiillt, falls | | < 2° gilt. Eine von F.W. SCHAFKE [12] hergeleitete Konvergenz- 
bedingung fiir die Potenzreihenentwicklung (0.3) ergibt die Forderung | ¢|<10. 
Durch genauere Ermittelung von 6 =f, laBt sich hier jedoch als Konvergenz- 
bedingung fiir das Iterationsverfahren auch |¢|<10 erreichen. Fiir Vektoren y 
ne 
C ist Ap =e| a und 


x7 


mit (7, Wo) =0, also Vektoren der Gestalt x=| 


0 
as e( 
|| 1X — (Ar x, Po) Poll | ¥? | Pol 


; te 
|Aiz|| bei vo=| 02. 


Die zweite der Ungleichungen (4.1) gilt daher fiir Elemente y der speziellen 
Form (4.2) p= B(A, p) — B(u,n) =— A, x (y= —7) fiir B5=B =| y@|, im Spezial- 
fall Yo=Yp also fiir B=0. Mit B=0 folgt aus (7.5) fiir Ayp=19, Po=Yo die Kon- 
vergenzbedingung | ¢|<10. 

Verwendet man B=|q| bei der Fehlerabschatzung (7.6) fiir die Aufgabe 
(10.1), so stimmen die Fehlerschranken mit den exacten Fehlerbetragen (Tabelle 4, 
Spalte 3) iiberein. Hieraus geht offenbar hervor, daB in § 4 wie auch in der 
vorangehenden Arbeit [1/6] vorsichtig abgeschatzt wurde. Praktisch wird man 
natiirlich mit B = 1+ ||q@—Yol| arbeiten. 


10.2. Beispiel 2 ist eine ahnliche Matrizenaufgabe 


(Mo+M)p=19, wvp=1 


; LOREaO Oi 
mit M,=(‘5 feist M, = «(‘ ‘| unde ei Of 
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Deny Werten y= 10, v)— (3 ,,benachbart ist die Lésung 


a 2* = e1(— 10+ /100 + e). 


A*(e)=10 + ez, g*(e) =| 


TAs — Ue w=(. erhalt man 
0 


Anti =10+ €2,,, m=. | (== OF AR eae) 
n 


wobel 


Znti=H(%,) ("= 0,1,2,...) mit f(2) = (1— 2) 


ist. Gilt 
eet) (4 (10.2) 


fiir die Lésung z* der Gleichung z=/(z), so konvergiert die Folge z, gegen 2*, 
wenn man mit einem 2) aus einer geniigend kleinen Umgebung von z* beginnt. 
(10.2) bedeutet 


é<10-/3, (10.3) 


Wahrend also die Potenzreihen von A*(e) und g*(e) den Konvergenzradius 
£22=10 besitzen, konvergieren die Naherungsfolgen der Iteration bei geeigneter 
Ausgangsnaherung gp noch unter der Bedingung (10.3). 

Diese Konvergenzbedingung kann man auch mit dem Satz aus § 6 erhalten. 
Es handelt sich hier um eine Aufgabe der speziellen in 7.2 untersuchten Art. 
Ebenso wie im Beispiel 1 la8t sich B=|qo| verwenden. Die GréBen a, b,c, B 
in (7.4) unterscheiden sich beliebig wenig von = |(A,@*, Po) |, 0, O bzw. | p*?| = 2%, 
wenn man Ay= My setzt und gp so wahilt, daB 


lv* — poll < 6 (10.4) 


mit geniigend kleinem 6 ist. Auf den nicht schwierigen Beweis sei hier verzichtet. 
Gilt also 


1—a2* > |(A, p*, Yo)| («=—), (10.5) 


so sind bei geniigend kleinem 6 in (10.4) auch die Konvergenzbedingungen (7.5) 
erfiillt. (10.5) fithrt aber gerade auf (10.3). 


§ 11. Beispiel einer Verzweigungsaufgabe 


Die §§ 8, 9 und 10 enthalten Anwendungen des Satzes aus § 6 auf Eigenwert- 
aufgaben in einem Hilbertschen Raum. Hier wird der Satz nun zur Lésung einer 
nichtlinearen Aufgabe in einem Banachschen Raum benutzt. 


11.1. Wir betrachten einen Stab, der an einem Ende eingespannt und am 
anderen Ende frei beweglich ist und welcher an dem freien Ende eine Last P 
tragt (Fig. 1). Die Auslenkungsfunktion g(x) (x = Bogenlange) geniigt (unter 
geeigneten Voraussetzungen) der Differentialgleichung 


— 9" =Apl1—(¢'? (1:44) 
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und den Randbedingungen di) = ote (14.2) 


wobei J bis auf einen Faktor die Last P bedeutet. 


Der Zusammenhang zwischen der Auslenkung 6= (1) und 4 wird (quali- 
tativ) durch Fig. 2 veranschaulicht (gestrichelte Kurven). Die entsprechenden 


i‘ of 


Fig. 1. Knicken eines Stabes Fig. 2. Zusammenhang zwischen Ausbiegung und Last 


Kurven fiir das linearisierte Problem mit der Differentialgleichung 


— 9" =Ay (11.3) 
sind in Fig. 2 punktiert gezeichnet. Diese vereinfachte Aufgabe liefert noch die 
richtigen Verzweigungsstellen (Knicklasten). 

Wir wollen hier die Aufgabe”° 
— P= hyp — Zh plyy, (11.4) 
PO) =97)=0, gil)=d>0 (11.5) 
bei gegebenem 6>0 lésen. Die Differentialgleichung entsteht aus (11.1), wenn 
man von der Taylorentwicklung der Wurzel nach (p’)? die ersten beiden Glieder 
berticksichtigt. Gesucht sind also A und (x) bei gegebenem 6=q/(1). Dabei 


gehen wir von der entsprechenden Aufgabe mit der Differentialgleichung (11.3) 
statt (11.4) als ungestértem Problem aus. Dieses besitzt 


os P= Y= 6-sin—= x 

a : 2 

als Lésung. Eine benachbarte Lésung der gegebenen Aufgabe (11.4), (44.5) soll 
berechnet werden. 


== 


11.2. Zunachst werde die gestellte Aufgabe funktionalanalytisch formuliert. 
Es sei © der Banachsche Raum der im Intervall [0, 1] stetigen Funktionen (x) 
mit der Norm?2 


loll =— Max | (x)| (11.6) 


00<x 


© Auf derartige Aufgaben wendete Péscut in [7] die Stérungsrechnung in etwas 
anderer Weise an. 

21 Man konnte ohne weiteres ||p|| = Max | p(*)| benutzen. Mit der Norm (11.6) 
lassen sich manche der auftretenden Gleichungen nur etwas bequemer schreiben. 
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und %& die Menge der Funktionen m€ © mit stetiger zweiter Ableitung, welche 
den Randbedingungen (11.2) geniigen. Mit den auf %& definierten Operatoren 


Ayp=—@", A,y=9' 
und 


B(A, y) =C(A, vy, 419) =— 349 (4,9)? GG) 
sowie dem auf © durch 


1 
Lp =~ (1) 


erklarten Funktional/ laBt sich die Aufgabe (11.4), (11.5) in der Form (1.7), 
(1.8) schreiben. 


Die Aufgabe (1.2) 


D) 


— 9" —*-~=(x), (0) =9'(1) =0 


ist fiir alle stetigen Funktionen w(x) lésbar, fiir welche 
1 
®y=0 mit Dy = { sin Sx p(x) dx 
0 


gilt, und die allgemeine Lésung hat die Form 


p(x) =cy(x) + Ry 


mit beliebiger Konstanten c und 


? sin x cos ZE fiir «sé 
TU 2 2 


Ry=fK(x6) p(6)dé bei K(x,8) = 


Doe p. 
*sin™ Ecos x fir <x. 
MA 2 2 


Es ist 


1 
LRy == (Ry)rn1 = 0 und Oy, =O f sin? xdz=— >0. 
0 


Ferner berechnet man mit partieller Integration unter Benutzung der Rand- 
bedingungen 


1 
72 


P(Ap— pol) p = f sin Fx gy 7 9) ax =0. 
0 


Die in 1.1 fiir den allgemeinen Fall genannten Voraussetzungen iiber Ay sind 
daher mit den oben definierten GréBen /, ® und R erfiillt. 


Man schatzt ab: 
1 ! 1 ~ vr 1 vr +: 
| Avgl|=-5 Max|o'|= 5 Max| Sg (6) a| <5 Max|¢ (6)| =||4ogl| fir pee, 


kann also in (1.10) Q,(s, 4)=¢ verwenden. Jedoch werden wir diese Funktion 
bei der Abschatzung nicht benutzen. Durch Ausmultiplizieren der auftretenden 
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Produkte C(...,...,...) stellt man fest, daB 


|C(A, p, Ary) — Clu p, Ar)| 
=|C((A—p) + (U— Mo) +o, (P—Y) + (Y— Yo) + Yo 
(A,p — Ay) + (Ary —A1 Po) + Ay Yo) — 
C((u— Lo) + Ho» (PY — Yo) + Yo» (Ayp—A, Ho) +41 ¥0)| 


<C(|A—p|+]4—ol +] Hol, 


g—yl+|y—Yol +|¥ol. 
|A.e—Aiy| +| Ary Ay Yo| +| 41 Pol) — 
, | Ary — 410] +] 41 Pol ) 
» O[Ile—vll+|l¥— voll +Ilvoll], 

6[||41¢— 41 || + ll41y — Ar voll + |] 41 Poll] ) — 
C(|e— Hol + lool, OE lly — voll + llyoll], OLl|41y — 41 ¥0l] +|] 41 Y0ll]) 


C(|u— Mol +] Hol, |Y—Yol +] Yo 
<C(|A—p|+|4—~Hol +] 40 


gilt?*. Die Beziehungen (1.11) sind fiir die Funktion (11.7) demnach bei w= wo 
erfiillt mit 


C= 4 (e+ Iual) 6 (ot + Ilvoll) [8(¢8 + |] Arvoll)?— 5 [aol (Allyl!) [|] Arvo ll. 


dhe 


4 


(11.8) 


Ee, ot, oh = (e+ Z)@+0 (e+ Z)-= 


11.3. Fiir Ap= Wo, Po= Yo ergeben sich A, und g, aus den (6.1) entsprechenden 
Gleichungen 


+ m2 ‘ nN / 3 1 12 
A, =o+ ty mit => | sin 5 x(—) 7) dsin 3 x|(dsin 2 a) dx, 
0 
i Hee Ne. fimerue Ae eam ee a Boole Gu mw og 3H, 
bees P= (Ay—Mo) dsin 5 x—— “7 Osin oa |(osin z x] | ren sin x, 
71 (0) = vi (1) = 0, pill) = 0 
zu 
Foes eens (x) = dsin 7 x — = 6(sin 2 « + sin 3% x 
eomee| La Pi a 2 256 | 2 ae 


11.4. Wir ermitteln die in der Fehlerabschatzung benétigten Konstanten. 
Esist 


mith (== + (V2 —1), 


2 In ahnlicher Weise kann man oft eine Majorante C ermitteln, wenn C ein Poly- 
nom seiner Variablen ist. 
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|A,Ry|| = =- y Max 


ie K(x, 6) p(6) dé) s 


ae K(x,8)|a8- |p| =callyll 
matonigg = (/2— 1), 


2, 
2.0 


fo 


See 4 
Dy, S2fsinZ xdx|ly||= Bolly] mit Po=h.=—, 


i 
f sin . x p(x) dx 


0 


ly — ay, Pv-val <5 Max(|woni +] 5 Join Zey@ ae|sin 2 2) 


S|lpl|-+2fsin Fédé-llyll=Allyll mit 6.=A.=14 4 


at 
| AL — Mol ee 
3 
lle: — Yol| = 8 Max |sin > x + sin =s x| = ae “3 &=o1, 
I|A1 (G1 — Yo)|l = =a 7 0 Max ‘cos 5 x J 3.cos om x| = a /2 = Or 


Zur Fehlerabschatzung sind die Gleichungen (6.3) mit diesen Konstanten 
a;, B;, 0, und 6} =0 sowie der Funktion C aus (11.8) zu lésen. Es wird t!—o! = 


Ps (x°—o?). Man kann also 7! eliminieren und erhalt dann leicht eine alge- 


0 
braische Gleichung 4. Grades fiir die GréBe a = se (x°—o}). 


Damit ergibt sich nach dem Satz aus §6 die Fehlerabschatzung fiir 6= ,4: 
| A*—A,| =| A*— wo— | =| A*— 2,467 40 — 0,007 61| =| A*— 2,475 01| S$ 0,00024, 


| p*(x) — gr ( 


*(x) — 0,099961 sin = x — 0,000039 sin > x| < 0,000008, 

fat OS: 

et ee 
| e*(x) — a(x) = |p*(x) — 0,19969 sin = x — 0,00034 sin 3™ 5 | S0,00028. 
Eine formale Potenzreihenentwicklung von a nach 6 liefert a = 2,234 + Glieder 


mit héheren Potenzen, also 


7° — of = 2,234 « 64 


-_ Glieder mit hdheren Potenzen. 
t! — 0; = 0,670- 64 


Die Potenzreihenentwicklung der Fehlerschranken fiir |A*—A,|, ||~*— || be- 
ginnt also mit der Potenz 6*, so daB — bei Konvergenz der Potenzreihen — das 
Verhaltnis der Fehlerschranken zur Anderung |A,—49| bzw. ||g,—q@ol| mit 46 
wie 6? gegen Null strebt. 
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[6] 


= 
ey 


[8] 
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Note on a paper ‘Further Remarks on the 
Stress-Deformation Relations for Isotropic Materials” 


A. PIPKIN & R.S. RIVLIN 


In § 114 of the above paper, the result is obtained that any polynomial scalar 
invariant under the group of orthogonal transformations of any number of 3 x3 
matrices is expressible as a polynomial in traces of products of the matrices. 
This result is in general valid only if all of the 3 x3 matrices are symmetric. 
More generally, if the matrices may be non-symmetric, the argument given in 
the paper leads to the result that any polynomial scalar invariant under the 
group of orthogonal transformations is expressible as a polynomial in traces of 
products formed from the matrices and their transposes. 

If one of the matrices is non-symmetric and the remaining matrices are sym- 
metric and the polynomial scalar invariant is linear in the elements of the non- 
symmetric matrix, then it may be expressed as a polynomial in traces of products 
formed from the matrices. This is the theorem which is used in §3 of a paper 
by A. E. GREEN & R. S. RIviin**. 


The results presented in this note were obtained in the course of research sponsored 
by the Office of Ordnance Research, U.S. Army, under Contract No. DA-19-020- 
ORD-3487. 


* J. Rational Mech. Anal. 4, 681 (1955) 
xx Arch. Rational Mech. Anal. 1, 1 (1957). 


Brown University, 
Providence, Rhode Island 
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Errata 


Correction to A. E. GREEN & R. S. RIVLIN, The Mechanics of Non-Linear 
Materials with Memory, Vol. 1, 1—21 (1957) 


Equation (3.12) does not follow strictly from the definition of F, since although 


Di; is an arbitrary tensor, B;; is not. We may, however, obtain equation (3.12) 


by the following argument. 
Since F may be considered either as a linear function of B,, or b;;, we have, 


with (3.10), 
(se Ox; 1 dis) a ay Ls 
ON MONE gt Kil) Tp ie OB ae bee 


oF Oe Sore ee O4; 1 »;) oF 
OX, OX, 7B Sek 4 ep 


Thus, in equation (2.28), we can without loss of generality take /,, to be given 


Baby. 3.19). 


(Received September 22, 1957) 


Erratum pour le mémoire de PHAM Mau Quay, Inductions électromagnétiques 
en relativité générale et principe de Fermat, Vol. 1, 54—80 (1957). 


Au 2éme théoréme du ICH Oe Sy, Lobes 
Dans un milieu isotrope de constantes diélectrique et magnétique e, u variables, les 
vayons électromagnétiques sont des géodésiques de longueur nulle de l’espace riemannien 
V, muni de la métrique 
a = 1 
aS = By pA xd xP = (cap = (1 — eal Matt) d x* d xB 


OU Lup est le tenseuy métrique fondamental et ug le vecteur vitesse unitaive d’univers 
définis en chaque point du milieu considéré. 


(Manuscrit regu le 7. octobre 1957) 
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